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Abstract: We present limit theorems for a sequence of Piecewise Determin- 
istic Markov Processes (PDMPs) taking values in a separable Hilbert space. 
This class of processes provides a rigorous framework for stochastic spatial 
models in which discrete random events are globally coupled with continuous 
space-dependent variables solving partial differential equations, e.g., stochastic 
hybrid models of excitable membranes. We derive a law of large numbers which 
establishes a connection to deterministic macroscopic models and a martingale 
central limit theorem which connects the stochastic fluctuations to diffusion 
processes. As a prerequisite we carry out a thorough discussion of Hilbert space 
valued martingales associated to the PDMPs. Furthermore, these limit theo- 
rems provide the basis for a general Langevin approximation to PDMPs, i.e., 
stochastic partial differential equations that are expected to be similar in their 
dynamics to PDMPs. We apply these results to compartmental-type models 
of spatially extended excitable membranes. Ultimately this yields a system 
of stochastic partial differential equations which models the internal noise of 
a biological excitable membrane based on a theoretical derivation from exact 
stochastic hybrid models. 
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1 Introduction 



In this study we present limit theorems for sequences of Piecewise Deterministic Markov Processes 
(PDMPs) with values in a separable Hilbert space. PDMPs are a particular class of cadlag, 
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strong Markov processes which combine continuous deterministic time evolution and discontinuous, 
instantaneous, random 'jump' events. We note that in view of applications this paper is ultimately 
motivated by the interest in the derivation of a justifiable Langevin approximation to spatio- 
temporal stochastic hybrid models of excitable membranes, e. g., neuronal membranes. This is 
accomplished by the limit theorems we present in the following. 

We start briefly introducing the general idea of our framework and the main results which are 
made precise in the subsequent sections. We consider a family of fully coupled, Hilbert space- 
valued PDMPs indexed by n € N. Here fully coupled means that the PDMPs which split into 
a continuously moving and a piecewise constant component are such that the jump rates of the 
processes depend on the state of the full system and the continuous dynamics depend on the state 
of the jump component. For the limit theorems we rely on two key assumptions. Firstly, jumps 
possess heights decreasing to zero for n — > oo but occur at an increasing frequency roughly inversely 
proportional to the jump size. We are therefore in the fluid limit setting, cf. [29j [30] . Secondly, 
we assume that for each n the continuous dynamics in between jumps depend on the piecewise 
constant component only via a finite set of (Hilbert space-valued) functions thereof, which we 
call coordinate functions. It is the sequence of coordinate functions coupled to the continuous 
component for which we derive limits. The first limit theorem we present is a weak law of large 
numbers for PDMPs in infinite-dimensional Hilbert spaces where the deterministic limit is given 
by a solution of an abstract evolution equation. Next we proceed to the presentation of a central 
limit theorem for the martingales associated with a PDMP. This central limit theorem gives the 
basis for an approximation of PDMPs by diffusion processes which are solutions of stochastic 
partial differential equations. Finally, we show how to represent the stochastic process arising as 
the limit in the central limit theorem as a solution of a stochastic partial differential equation 
(SPDE) which then yields a Langevin approximation for PDMPs by a system of SPDEs. The new 
results presented extend previous results for PDMPs and pure jump processes in Euclidean space 
[501 IT51 155] . The difficulties in extending the fluid limit theorems in 133 EH] to processes taking 
values in infinite-dimensional Hilbert spaces lie, on the one hand, in the appropriate treatment of 
Hilbert space-valued martingales. These arise by splitting a PDMP, being a semi-martingale, into 
a sum of a part with finite variation and a local martingale. As these considerations are essential 
we have devoted a full section, Section E2 to the discussion of the martingales. On the other 
hand, the more intricate existence theory of solutions to abstract evolution equations compared 
to solutions of ordinary differential equations in Euclidean space demands for additional technical 
rigour. 

We apply our theoretical findings to spatially extended hybrid models of excitable membranes. A 
first hybrid formulation of one such model in the context of neuroscience was presented in [5] and 
reformulated and extended as examples for PDMPs taking values in infinite-dimensional Hilbert 
spaces in [T3|. For example, the Hodgkin- Huxley model is a deterministic, macroscopic model for 
the coupled evolution of the neuronal membrane potential and the averaged gating dynamics of ion 
channels [21] . More realistically, the membrane potential, which is the macroscopically observed 
variable of interest, arises from the stochastic dynamics of finitely many ion channels. Thus the 
application of our limit theorems shows that the Hodgkin-Huxley is obtained as the limit of a 
sequence of stochastic microscopic models taking the form of Hilbert space valued PDMPs in the 
sense of a law of large numbers. Conceptually, here the fluid limit corresponds to increasing the 
number of ion channels while simultaneously decreasing the individual influence of an individual 
channel on the total current. The martingale central limit theorem can then be used to define the 
Langevin approximation providing a relatively simple stochastic version of the Hodgkin-Huxley 
model incorporating internal fluctuations. 

Concluding this introduction, we comment on related work to fluid limits in the infinite-dimensional 
setting. Averaging for PDMPs in infinite-dimension, in particular for the neuron model introduced 
in [5], wherein also a law of large numbers was considered, has been recently considered in [20] . 
For a model of linear chemical reactions by jump Markov processes a law of large numbers [4] 
and a central limit theorem [35] have been proven based on the original work of [5^1 130] for finite- 
dimensional jump-processes. In these cases the deterministic limit is a reaction-diffusion partial 
differential equation and the central limit theorem yields diffusion processes given by stochastic 
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partial differential equations. Limit theorems for variations of this model have been investigated 
in two series of studies, cf. [26j [27l [28] and [7j [8j [9l [TUJ [11]. A central difference between spa- 
tial models of excitable media to models of chemical reactions is that the latter exhibit diffusive 
motion of the reactants (~ channels) which is absent in the former. Additionally excitable me- 
dia equations exhibit non-local interaction of channels as their dynamics are coupled globally via 
the membrane potential. The limit theorems we establish have to account for these differences. 
Further, there is also a difference on the technical side. The technique employed in [25] and in 
all subsequent publications cited above is based on the semigroup approach to stochastic / de- 
terministic evolution equations. In contrast, we pursue in the present paper the approach of a 
weak formulation. At large, the weak formulation of evolution problems allows to consider more 
general equations as when dealing with mild, strong or classical solutions, cf. a discussion of this 
aspect in [321 Chap. 23.1]. Finally, we also mention a central limit theorem for Hilbert- valued 
martingales in [34] and a diffusion approximation of SPDEs on nuclear spaces driven by Poisson 
random measures in [23) . The methods of proof we employ for the theoretical results in this study 
are motivated by the two last references, but differ as the classes of stochastic processes considered 
therein and in the present manuscript are different. 

The remainder of the paper is organised as follows. We first briefly define PDMPs in Section [2] 
and precisely state the structure for a sequence of such PDMPs to allow for a limit. Then we 
discuss in detail the associated martingale process in Section [3] Limit theorems and the diffusion 
approximation are presented in Scctions[3]and[5] We have deferred the proofs of the main results to 
Section^ Next in Section[7]wc discuss applications of these limit theorems to compartmental-typc 
models of excitable membranes where the proofs of the conditions are deferred to Appcndix[B] The 
paper is concluded in Section[8]with a brief discussion and an outlook on further developments and 
applications. Finally, the Appendix [X] of the paper contains the proof of the technical Theorem 
I3.1l that guarantees the square-integrability of the associated Hilbert space valued martingales and 
establishes an appropriate Ito-isometry. 

2 Piecewise Deterministic Markov Processes 

In the first subsection we briefly define PDMPs and, in particular, discuss the specific subclass of 
PDMPs for which we present limit theorems in this study. For a general discussion of PDMPs we 
refer to the monographs [THl [22] and, specifically, for Hilbert space valued PDMPs associated to 
solutions of partial differential equations we refer to [131 HO] • In the second subsection we present 
the sequence of PDMPs for which the limits are analyzed in this study. Finally, a notational 
remark: in this paper pairings (-,-) and (-,-) denote the inner product or the duality pairing, 
respectively, with respect to a certain Hilbert space which is usually indicated with a subindcx. 
Further, * is used to denote dual spaces. 

2.1 PDMPs on Hilbert spaces 

Let (fi, J", {J-t)t>o,^) denote a filtered probability space satisfying the usual conditions, X C H C 
X* be an evolution triple of separable real Hilbert spaces and K be a countable set of isolated 
states. The product HxK serves as the state space for a PDMP. Then, a PDMP is a cadlag strong 
Markov process X t (co) = (U t (uj), Ot(w)) £ H x K for all t > which consists of two components. 
The first, Ut, takes values in H, possesses continuous sample paths and is denoted the continuous 
component of the PDMP. The second, 0t, taking values in K and possessing right-continuous, 
piecewise constant sample paths, we call its jump component. We say a PDMP is regular if the 
number of jumps of 9< is a.s. finite in every finite time interval [0,T]. In this study PDMPs are 
always regular. 

We next state the mechanisms which govern the time evolution of the paths of such a PDMP. 
Firstly, there exist for each 8 £ K an abstract evolution equation 

u = A{9) u + B($, u) (2.1) 
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where A(9) : X — > X* is a linear and B(9, ■) : X — > X* a (possibly nonlinear) operator. We assume 
that the family of abstract evolution equations (|2.1j) is well-posed, i.e., given any 9 £ K and any 
initial condition u £ H there exists a unique global weak solution <j)(-,(u,0)) £ L 2 ((0,T), X) n 
i/ 1 ((0,T),X*) depending continuously on the initial condition. Note, that the regularity implies 
4>(-,(u,8)) £ C([0,T],H), cf. [351 Chap. 11]. Then the trajectory of the continuous component 
Ut follows in between jumps of the jump component 0t the weak solution to (|2.ip corresponding 
to the parameter 9 given by the current state of the jump component. That is, for T k , k £ N, 
denoting the jump times of the PDMP we have that 

U t = <f>(t - T k , (U Tk ,Q Tk )) Vt £ [Tfc, Tk+l) ■ 

Secondly, describing the stochastic transition dynamics of the jump component 0t there exist 
measurable transition rates A : H x K R + that define the distributions of the random jump 
time of t in the sense that for all 9 £ K 

P[0 t+S = t , O<s<At|0 t ] = vq>(- J A(U t+s , t ) ds) . (2.2) 

In view of (|2.2|) we assume that A is intcgrablc along the solutions of (|2.1|) on any finite time 
interval, i.e., 

/ A((j){t,(u,9)),9)dt < oo VT<oo 
Jo 

for all 9 £ K and all initial conditions u £ H. but diverging as T — > oo. We note that in 
applications we usually find that the transition rate A is bounded which implies the regularity of 
the PDMP. Finally, there exists a Markov kernel fj, on H x K into K that gives the distribution 
of the post jump value, i.e., 

P[0 4 = e|0t^0 4 _] =/x((U t ,0 t _) ! {e}) VCeif. (2.3) 

The elements of the quadruple (A, B, A, /i) arc called the characteristics of the process and under 
the above conditions define a regular PDMP uniquely (up to versions). Furthermore, under 
these conditions the following result characterising the extended generator of PDMPs is proven in 

muni. 

Theorem 2.1. A function f : H x K — > R is in the domain of the extended generator of 
a PDMP if the mapping t i— > f(Ut,Qt) is absolutely continuous almost surely and the map- 
ping (£, s,u>) i— > f(U s —,£) — /([/ s _(w), s _(w)) is integrable with respect to the random measure 
A(C/ s _,0 s _) M ((f/ s _,0 s _),de)ds. 

Moreover, if in addition f is continuously Frechet- differ entiable with respect to its first argument 
such that the Riesz Representation f u £ H of the Frechet derivative satisfies f u {u,9) £ X for 
u £ X and is a locally bounded composition operator in L 2 ((Q, T), A")Jf| then the extended generator 
Af is given by 

Af(u, 9) - (A{9) u + B(9, u), f u (u, 9)) x + A(tt, 9) (/(«, £) - f(u, 9)) fx((u, 9), d£) . (2.4) 
2.2 An appropriate sequence of PDMPs 

Let E denote another separable real Hilbert space. Further, for a certain m £ N (its significance 
is explained in the next paragraph) we denote by % = X i=1 H, £ = X i=1 E the direct products 
of the Hilbert spaces H and E which are Hilbert spaces themselves. Finally we setf* = X-LiE* 
which is the dual space to £. 

1 Note that the Frechet derivative at a point u S H is a linear, bounded functional on H and thus an element of 
the dual H * . 

2 An example of such a function / is (u, 8) i— > | u, 1 1 ^ in which case fu(u, 9) = 2u. 
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We now define the structure of the sequence of processes for which we derive the limit theorems. For 
all n £ N let (0™, F n 1 (J r ")t>o,P n ) be a a filtered probability space satisfying the usual conditions 
and the processes (X t n ) t >o = (£7™, 0™)t>o defined thereon are regular PDMPs taking values in 
H x K n with path properties as defined in Section |2~T1 Correspondingly, the characteristics of the 
PDMPs are given by (A n , B n , A™, fj, n ). Note that the state space K n for the piecewise constant 
component changes with varying index n whereas the state space H for the continuous component 
remains fixed. Therefore, in order for such a sequence of processes to allow for a limit we need to 
impose a special structure on the characteristics referring to the continuous component. To this 
end we assume there exists an m <G N, introduced above, such that for each PDMP (U™, 0™)t>o 
there exists a family of measurable coordinate functions zf : K n — > E, i = 1, . . . , m, such that the 
characteristics A n (0) , B n (6) depend on the piecewise constant component and on the index n only 
via the £ -valued coordinate process z n (6) = (z™(6) , . . . , z^(#)). That is, there exist measurable 
operators A, B : £ x X — > X* such that for all n € N, all u <G H and all e K„ 

A n {6)u = A(z n (6))u, B n (6,u) = B(z n (6),u). (2.5) 

The coordinates z n can be interpreted as a 'sufficient statistic' of the piecewise constant component 
for the evolution of the continuous component. In statistics a sufficient statistic for a quantity of 
interest is a function of the observations that is sufficient to estimate this particular quantity. For 
example, the sample average of independently and identically distributed real random variables is 
a sufficient statistic for the mean of their distribution. In the present setting, this means that the 
coordinate functions contain all information about the vector 9 that is needed to determine the 
continuous dynamics in between jumps. Further, the essence of the subsequent limit theorems is 
that the sequence of coordinate processes on the space £ allows for a limit under certain conditions. 
Typically, in applications one is interested in the dynamics of the continuous components only, thus 
a restriction of the attention to the coordinate functions is well justified. As £ is a (vector-valued) 
Hilbert space itself, no generality would be lost if instead of the family of coordinate functions 
we assumed the existence of Hilbert space-valued functions z n taking values in the same Hilbert 
space for each n. However, we decided to use this more detailed notation since in examples one 
usually encounters that it is a set of coordinate functions that encodes the information necessary 
for defining the dynamics of the continuous component. 

In order to illustrate this set-up let us briefly discuss the Hodgkin-Huxlcy model as an example 
of the general excitable membrane model considered in Section [7] Here the sequence of abstract 
evolution equations (|2.1[) arises from parabolic partial differential equations modelling the space- 
time evolution of the membrane potential of the form 

u(t, x) = Au(t, x) + 9iPi(t,x)(Ei-u{t,x)), (>0,i£flcK' i (2.6) 

i=Na,K,L 

with constants g i > and Ei £ K, cf. (|7.6[) . The indices refer to electrical currents due to 
the movement of charged Sodium (Na) and Potassium (K) ions across the membrane and ohmic 
leakage (L) current mainly due to Chloride ions [2U- In hybrid versions of the Hodgkin-Huxley 
system the conductances pi (t, x) depend on the finite number of open ion channels distributed 
in the membrane which increases with n. Each individual channel is modelled stochastically 
opening or closing at random times with dynamics depending on u, cf. Section 17.21 for more 
details. In the case of constant potential u(t, x) = u each channel were a continuous time Markov 
chain. The collection of channel states at any time instant t defines the discrete component O™. 
Finally, the coordinate functions z n relate channels in a specific state to their location in the 
physical space D, cf. their definition in (|7.5[) . They map the channel configurations into piecewise 
constant space-time functions stating the local density of channels in the particular states, thus 
p™(t,x) := z™(9") <G L 2 (D) = E. Hence, equipped with suitable boundary conditions equation 
(|2.6[) is an abstract evolution equation of the type (|2.ip where the Hilbert spaces H, X and E are 
spaces of real functions onDc M. d . 
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3 The associated martingale process 



For the limit theorems we derive in this paper, the main estimation procedures concern certain 
martingales associated with the PDMP. As these are of such central importance we discuss them 
in this separate section. The principle aim is, on the one hand, to derive conditions that imply the 
convergence in probability of the associated martingales as needed for the law of large numbers 
(cf. condition (|4.5[) in Theorem 14. ip and, on the other hand, we present some necessary structure 
for the central limit theorems. Therefore we define for all j = 1, . . . ,m the E- valued stochastic 
process M™ by 



M?(t) := 4(e?) - *;(es) - / [A n {- , *;(•)>*] (u?, e; 1 ) d*, (3.1) 

Jo 

where the integrand in the right hand side is given by 

[A n (-,z^(-)) E ](ui\Q^ = A"(E/?,e?) f -*;(e?))^ B {(u?, e?),df) 



K„ 

= k'\u:\ &:) (*;(o - 4( e ")) ^((^ e ")> ui) ■ ( 3 - 2 ) 

Hence the integrand is a countable convex combination of elements in E with time-dependent 
coefficients and in between jumps it depends continuously on s. Anticipating condition (|3.4p 
below, which we generally assume to hold, we find that the integral in the right hand side of (|3.1[) 
almost surely exists in the sense of Bochner. For a brief discussion of the Bochner integral we refer 
to [37J App. A]. For an application of a functional 4> <G E* to (|3.ip we obtain 

(0, Mf{t)) E = (0, z?{Q n t )) E - (0, b - [ [A n (<P, *?(•)>*] 0J) ds, (3.3) 

Jo 

where the integrand is 

H n <& *;(•)>*] (t?,e?) = A"(c/ S ",e; 1 )/ - (0,4 l (e^)) £/ i"((c/ s ",e:),de). 

Thus the integral has the form of the extended generator, cf. Theorem 12. 11 applied to the mapping 
(w, 0) i-> (</>, Zj-(8)) e- This already suggests that the processes Q3.3P are martingales under suitable 
boundedness conditions. In fact we are able to establish that the processes M™ are E'-valued cadlag 
martingales. We refer to |14[|37j for a brief discussion of martingales in infinite-dimensional spaces. 
The easiest way to validate the martingale property is due to the following result (37J Sec. 2.3]: 
If E"||M™(4)||e < oo for all t G [0,T], the Hilbert space-valued martingale property holds if and 
only if (0, M"(t)) E is a real- valued martingale for all S E* . The following theorem gives a 
condition that guarantees that the processes (|3.ip are square-integrable martingales and satisfy 
an Ito-isometry. The proof is rather technical and thus we have deferred it to the Appendix [A] 

Theorem 3.1. Let n G N be fixed and assume that for all t > it holds that 

e /"V(E?,e?) / - ^(e?)|||A* n ((c^,e?),d^)l ds < oo. (3.4) 

TTien £/ie process M" is a square-integrable martingale and satisfies the Ito-isometry 

E"||M;(t)|||= /V[A"([/;,e?) / u*?® -^(e?)|||. M n (0C, e?), de)l ds. (3.5) 

We continue the investigation of the processes M™ as Hilbert space valued martingales. From 
now on we always assume that assumption (|3.4p holds. Note that the finiteness of the second 
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moments of the jump sizes is a standard condition in related fluid limit theorems [5^1 [35J l3~i] . 
We introduce a concept akin to the quadratic covariance operator in Euclidean finite dimensional 
spaces. This concept is important for the central limit theorems in, on the one hand, establishing 
weak convergence, and, on the other hand, characterising the limit. For further reference we refer 
to [33]. 

Definition 3.1. For the square-integrable, _E-valucd, cadlag martingale M™ we denote by (<C 
M"3> t ) t >o its predictable quadratic variation process, i.e., the unique (up to indistinguishability), 
predictable L\(E* , _E)-valuccjf| process which satisfies that for all (f>, ip € E* the real- valued process 

t ^ (0, M?(t)) E (il>, Mf{t)) E - (0, «A/J l »^) B (3.6) 

is a local martingale. 

The aim now is to obtain an explicit formula for the quadratic variation process of the individual 
martingales M™ as well as of the vector- valued process M n of all martingales M™, i.e., the £ -valued 
process 

t^M n (t) = (M?(t),...,M£(t)). 
To this end we define for alH, j = 1, . . . , m operators G™- G L(E* , E) by 

1>» G^(u,0 n )i/j := (3.7) 

Clearly, these are linear, bounded operators mapping E* — > E and depend measurably on (u, 9 n ) G 
H x K n . For i — j each operator is non-negative, i.e., (0, G!?.(u, > for all <f> G E*, and 

symmetric, i.e., (?/>, GJj(u, O n )<f)}E = (4>iGJj(u,8 n )ip)E for all (j), if) G -E*. Let (<£fc)fceN denote an 
orthonormal basis in E* . We find due to the Riesz Representation Theorem and Parseval's identity 
that the trace of the operators Gjj satisfies 

TrG^(u,n = A n (u,n / x)(^*.«i(o-«"(n>B)% n ((«»n,df) 

= A>,r) / ii^o-^cniil^CKn,^). (3.8) 

For arbitrary i,j the trace is bounded in terms of p.8[) as it follows from Young's inequality that 
TrG£(u,0") < ±TrG?t(u,6 n ) + §Tr G!?(u, 0"). 

Let $ = (0i, ... , m ) and ^ = ■ ■ ■ , V'm) be elements of £* . Summing over all operators (13.71) 
applied to the components of 4>, \1/ as indicated by the indices, i.e., 

m 

(<f>,G n (u,6 n )*) £ := £ {(j>i,G^{u, O n )ipj) E , (3.9) 

we obtain a linear, bounded operator G n (u,O n ) mapping £* to £. This operator is symmetric as 
the operators G™- satisfy (</>, G^(u, Q n )ip) e = (Vs G^{u, n )(j>) e for all i,j. Moreover, the operator 
G n (u, 9 n ) is non-negative as it holds that 

*?(0 " < l (0) B ) M™ ((«, 0") , dO . 

Finally, the operator G n (u, 9 n ) is of trace class if the operators Gjj, j = 1, . . . , m, are of trace 
class and the trace satisfies 

Ti-G n (u,6 n ) = A n (u,9 n ) [ \\z n (£) - z n (8 n ))\\ 2 £ ti n ((u,6 n ),d() . (3.10) 
We next prove that the operators (|3.8[) give the quadratic variations of the martingales (|3.1|) . 

3 Li(E", E) denotes the space of trace class operators from the Hilbert space E* into E. 
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Proposition 3.1. The quadratic variation of the martingale M™ satisfies for all t > 

<MJ l >t = f G™ ([/™,9™)ds. (3.11) 
Ja 

Remark 3.1. It is an immediate consequence of Proposition 13.11 that the quadratic variation of 
the £ -valued martingale M n is given analogously to p. lip by integrating the operator G n . 

Proof. First of all note that due to the characterisation of the trace Q3.8[) and condition (|3.4[) it 
holds that the process in the right hand side of (|3.11[) takes values in Li(E*,E) almost surely. 
Further, it holds that <M™» t satisfies for all 0, tp e E that 

(0,«Af;»^) £ = 

as this right hand side is, due to [HI Prop. 4.6.2], the unique real- valued process such that 
(</>,M?(t)) E {ijj,Mf{t))E-{^<- M j">t^)E is a local martingale. Here (<j), M"(t)) E and {^,Mf{t)) E 
are understood as real-valued stochastic integrals with respect to the associated martingale mea- 
sure of the PDMP. Thus we infer that for all (j), ij) € E it holds 

(0,«A/;» t V) B = f {<j>,Gl{U:,@ n s )i>) E ds. 
Jo 

Finally, the linearity of the Bochner integral (note that L\(E*,E) is a Banach space) implies 

(EU). □ 

A further second property of the quadratic variation is that the process 

t h-> — Tr <M">( 

is a local martingale. We note that the trace process t t— > Tr <CM™;^ t is the unique, predictable 
increasing process exhibiting this property. Using the characterisation (|3.1ip of the quadratic 
variation we thus obtain that the process 

t (->• ||M"(i)||! — Tr Gjj(U™, 8") ds^j = \\M?(t)\\% - j TrG^([/ s ™, 6?) ds (3.12) 

is a local martingale vanishing almost surely at t = and analogously in the case of the f-valued 
martingale M n . 

We are now in a position to state a lemma which establishes the convergence in probability (|3.14p 
of the processes (M")t>o necessary for the law of large numbers, cf. condition (14.51) in Theorem 
IP 



Lemma 3.1. Assume that for all T > 



lim E 

n— >oo 



A n (u:\e:) \\zno-z?(®:)\\ 2 E v n {(u:\e:),dt) d S = o. (3.13) 

1 JK n J 

Then the process (|3.12p is a martingale and for all T, e > 0, it ZioWs i/iai 

lim P"[sup te[0)T] \\M?(t)\\ E > e] = 0. (3.14) 
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Proof. As the process M" is an E'-valued cadlag martingale, it holds that ||M™|||; is a cadlag 
submartingale. Thus an application of Markov's and Doob's inequalities yield the estimates 

su Pte[0)T] ||M;(t)||i> e ] < - e E n [m Pte[0tT] ||m;(*)|||] < ^e"||m;(t)|||. 

Now, the Ito-isometry p. 51) and condition (|3.13[) imply the convergence in probability (|3.14[) . It 
remains to show that the process (|3.12[) is a martingale. A sufheient condition, see, e.g., [HI 
Prop. B.0.13], is that for all T > it holds 



su Pte[o,T] 



\M?(t)\\ 



TrG"(C/r,e:)d S 



< oo 



(3.15) 



Estimating the term inside the expectation we obtain 



su Pte[o,r] 



mm 



TrG"(C/ s ",e:)d S 



< su Pte[0iT] ||M«(t)|||+sup te[0iT] / A n (u,e n ) / \\z?(0-z?(Q n )\\ E (* n ((u,e n ),dOds 



K„ 



The expectation of the first supremum term in the right hand side is bounded due to Doob's 
inequality and the square-integrability of the martingale. The term inside the second supremum 
is increasing, thus its expectation is finite due to condition (|3.13l) . 

□ 



4 A weak law of large numbers 

In order to propose a deterministic limit of the sequence of PDMPs we consider functions Fj : 
£ X H — > E, j = 1, . . . , m. In combination with the operators A, B these functions are used to 
define a coupled system of deterministic abstract evolution equations 

ii = A(p) u + B(p.u), 

' ' (4.1) 

Pj = F j(P> u ), j = l,...,m. 

We assume that to suitable initial condition (uq,pq) <G H x £ there exists a unique weak solution 
(u(t),p(t)) t >o in C(M. + ,H x £) of (|4.1[) . Additionally, we assume that for all i = 1, . . . , m the 
components pi satisfy 

{4>,Pi{t))E = {4>,pM)E+ I (^FMs),u(s))) E ds vte [o,T], 4>&e*. (4.2) 

Jo 

That is, the components pj satisfy the equation (|4.1[) in the sense of an Hilbert space valued integral 
equation. We note that in application one usually encounters deterministic limit systems that 
possess strong or classical solutions and hence the current weak framework is satisfied. Finally, 
we assume that the operators A 7 B and Fj, j = l,...,m, satisfy Lipschitz-type conditions on 
L 2 ((0, T), £ x X) in the sense that for every T > there exist constants L\ and L2 such that for 
all u,ve L 2 ((0,T),A) and all p, q e L 2 ((0,T),£) it holds that 

T 

(A(q) v — A(p) u, v — u)x + (B{li v) — B(p, m), v — u)x dt 



and 



f 

< Li / \\v-u\\ 2 H + y>i-ft||^di. (4.3) 
J o i= i 

fT 2 r T m 

(j o \\Fi(q>v)-FjM\\Edt) <L 2 J^ Wv-utu + J^U-piWldt, 



(4.4) 



where we have omitted the arguments t of the functions u,v,p and q. 
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Remark 4.1. In the proof of the law of large numbers, see Section IP! these Lipschitz conditions 
are applied such that one pairing (v, q) refers to a path segment of the continuous component of 
a PDMP and the coordinate process and the second (u,p) to the deterministic limit functions. 
Thus for the applications of (|4.3|) and (|4.4[) in the proof it is sufficient that these hold only for 
pairings (v, q) out of a set containing almost all paths of the sequence of PDMPs and (u,p) being 
the deterministic limit, i.e., one (!) distinguished pairing. This restriction of (|4.3p and (|4.4j) to be 
satisfied only for particular pairings (v,q) and {u,p) out of the whole path space has a decisive 
advantage: In order to establish these conditions we are able to incorporate additional qualitative 
results on the trajectories of the PDMPs and the deterministic limit and the constants L\, L2 
may depend on (u,p). For example, in the application to excitable membrane models such an 
additional qualitative is that the components corresponding to u,v,p,q are pointwise bounded. 

We now present a weak law of large numbers in Theorem 14.11 below. The proof of the theorem 
follows the lines of previously published limit theorems considering processes in finite dimensions 
[23 EEI • The main difficulties arising in infinite-dimensional phase space concerns the bounds on 
the martingale part, cf. condition (CI), which is rarely a problem in finite dimensions. However, 
using the appropriate martingale theory in Hilbcrt spaces these can be kept to a minimum. Then 
the difficulties arc mainly of a technical nature as martingale theory in connection with PDMPs 
in infinite-dimensional spaces gets more involved and is not covered by previous results in [22]. 
We have established the necessary theory in the preceding Section [3] and addressed the question of 
the convergence of the martingale part (CI) within this framework. Most importantly, in Lemma 
13.11 we have proven a sufficient condition for (CI) to be satisfied. In particular, this sufficient 
condition (|3.13[) is a natural extension of the condition employed in finite dimensions, cf. [291 135] . 
A different approach to establishing condition (CI) which avoids using martingale theory in Hilbert 
spaces is exemplified in the law of large numbers proved in [5] ■ In infinite-dimensional space this 
approach encounters the problem of simultaneously controlling countably many real martingales 
compared to only finitely many in the case of its finite-dimensional counterpart. This problem 
can be overcome with an intricate compactness argument which relies on the assumption that 
the dual space E* is compactly embedded in some additional normed space and all estimates - 
especially an estimate which also implies condition ()3.13[) - have to be derived in the norm of this 
additional space. Furthermore, the condition, that all martingales (((j), M™(t)),g) t >o, j = 1, ■ ■ ■ ,m 
and cj) G E* , possess almost surely uniformly bounded paths, has to be introduced. We are of the 
opinion that our approach is more elegant, but, more importantly, it avoids the introduction of 
additional conditions. 

Finally, consistently with the notation in Section [3] we use in the subsequent theorem and its proof 
the notation [A n {■ , {■)) E ] as defined in $£2$. Then, for given (u,8 n ) G H x K n functionals 
[A n ( ■ , z] l {-)) E ] (it, n ) on E* are defined by the mappings (j) i-> [A n ((j), z™(-)) E ] (u, 8 n ). As usual 
we identify the bidual E** with E and thus [A n (-, z™(-)) e] (u, n ) G E. 
Theorem 4.1. We assume that the following conditions hold: 
(CI) For all j = 1, . . . ,m it holds that for fixed T, e > 



where we have omitted the argument t of the functions u and 6. 

(C3) The initial conditions (Uq,Qq) of the sequence of PDMPs converge in probability to the 
initial conditions of the deterministic limit in the sense that for all e > 



n 



hm P"[sup te[0iT] ||M;(t)|| B >e] =0. 



(4.5) 



(C2) 



The functions Fj , j — 1, . . . ,m, satisfy for all e > that 




m 



71 



lim P" \\UZ-u \\ H +Y / \\z?(O 1 3)-p i (0)\\E>e =0. 



i=l 
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Then, for every e > and every fixed T > it holds that 



Urn 



S u Pte[0 , T] (||c/ t " - u(t)\\ 2 H + E K( @ t) -ft Will) > 



(4.7) 



Remark 4.2. The result (|4. T[) implies convergence in probability of the processes (U™, z"(0™)) 4 >o 
to the deterministic function (u(t) , p(t)) te [ x] in the Hilbert space L 2 ((0, T),H x £). If the differ- 
ences of the components are almost surely bounded independent of n the convergence even holds in 
the mean, cf. the application of the law of large numbers in Theorem 1 7. II Further, the conditions 
(C1)-(C3) are generalisations from Euclidean space to infinite-dimensional Hilbert spaces of those 
employed in the corresponding theorems for PDMPs in Euclidean space [35) and, in particular, of 
the original formulation in case of pure jump processes in Euclidean space |29| . In these cases the 
conditions above reduce to the corresponding assumptions. 



5 The central limit theorem and the Langevin approximation 

We proceed to the presentation of the central limit theorem for associated martingales (M t n ) t >o 
defined in (|3.1|) . The central limit theorem provides the theoretical basis for an approximation 
of spatio-temporal PDMPs by Hilbert-space valued diffusion processes where the latter can be 
represented by solutions of stochastic partial differential equations, see Section [5.21 
Proving central limit theorems usually involves two tasks: On the one hand, one has to show 
the existence of a limit and, on the other hand, one has to provide a characterisation of the 
limit as a certain stochastic process. The former is equivalent to the problem of tightness of 
the stochastic processes. In the case of martingales sufficient conditions for tightness depending 
on the quadratic variation process are stated in |34j . In order to characterise the limit there 
exist different approaches, showing either that the limit solves a given (local) martingale problem 
which is known to have a unique solution (cf. |23[ I34j ) or proving weak convergence of the finite 
dimensional distributions (cf. |25l I35j ). We present two central limit theorems, Theorems 15. II and 
15.21 employing the two methods, respectively, however, to avoid repetition we state only the proof 
of the first in the present study and refer to the PhD thesis of one of the authors [ID] for the proof 
of the second. The two theorems differ in a technical assumption which in each case arises in 
addition to the central condition of the convergence of the quadratic variations. We believe that 
for applications of the limit theorems it is advantageous to know both versions of the martingale 
central limit theorem, as it is easily conceivable that only one of these technical assumptions is 
satisfied. Hence the theorems are applicable in different situations. 

Finally, we emphasise that in the following the space £ need not necessarily be the same space for 
which the law of large numbers is satisfied. However, clearly, the space £ in the present section 
contains the space in the law of large numbers as subspace. In applications, usually, the law of 
large numbers holds in a space with a stronger norm, for example, for the excitable membrane 
model considered in Section [7] the law of large numbers holds in L 2 (D) whereas the central limit 
theorem holds in the space This is a major difference to the corresponding results in 

finite-dimensional space where both limit theorems hold in the same spaced 



5.1 A martingale central limit theorem 

In this section we present central limit theorems for the scaled £-valued martingales {yfa n M™) t >o 
associated with a sequence of PDMPs where a n <E R+, n € N, is a suitable rescaling sequence such 
that limn-j.oo a n = oo. Clearly, the rescaling is necessary in order to be able to obtain a limit 
different from the trivial limit as (j4.5[) implies that (M™) t >o converges to zero in distribution. We 
note that the sequence a n can also be interpreted as characterising the speed of convergence of 
the martingales (M™) t >o. 

4 Here and everywhere else H~ 2s (D) is the dual space to the Sobolev space H 2s (D) where Del' 1 and s > d/2. 
5 Note also that in finite-dimensional spaces all norms, and hence also all norms on subspaces, are equivalent 
which does not hold in the case of an infinite-dimensional Hilbert space. 
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In the following let t i-> G(u(t),p(t)) G L {£*,£) be a Bochner-integrable operator-valued map 
such that each G(u(t),p(t)) is a symmetric, positive trace class operator. Particularly this implies 
for all $ G £ * and all t > 0, that it holds that 



($,G(u(s),p(s))$) £ ds< 



(5.1) 



Here (u(t) , p(t))t>o is the deterministic limit obtained in Theorem 14.11 and the use of this notation 
for the - at this point - arbitrary time-dependent operator G only illustrates that in applications 
the time-dependence is due to a dependence on the deterministic limit system. These operator- 
valued functions are used to define a unique centred diffusion process on £, i.e., an £ -valued 
centred Gaussian process with independent increments, continuous sample paths. Such a process 
is uniquely defined by its covariance operator and due to a theorem of Ito stated in |25j every 
family of trace class operators C*(t) G L\(£,£) which are increasing and continuous in t define a 
centred diffusion process. In the present situation we define C* in the following way. We denote 
by t : £ — > £* the canonical identification of a Hilbcrt space with its dual, hence we can define for 
x,y e£, 



,G*{t)y) t 



(i(x),G(u(s),p(s))t(y)) s ds 



which is continuous and increasing for all x £ £ and C* (t) is a trace class operator on £. Moreover, 
for operators C(t) G Li(£*,£), defined by 



($,(G(u( 8 ),p(sp)>, 



(5.2) 



there is obviously a one-to-one relationship between C* and C. Hence, we may say that also the 
latter defines a diffusion process on the space £ . 

We proceed to the statement of the central limit theorem. The proof of the theorem employs a 
characterisation of the limit via the local martingale problem. The essential condition character- 
ising the limit is the convergence of the quadratic variation processes (|5.6p . The second condition 
(|5.7[) is a technical condition on the jump heights which arises due to the method of proof and is 
usually satisfied in applications. The remaining conditions are such that (Dl) guarantees tightness 
of the sequence of processes and in combination with (D2) that any limit is a continuous stochastic 
process. The proof of the following theorem is deferred to Section [ 

Theorem 5.1. We assume that the following conditions hold: 

(Dl) For allt>0 it holds that 



sup a n E n 

n£N 



I r 



A»([/ S ",e; 1 ) / \\z n (o - z"(e?) Hi /i™ ((c/;, e?),de) d s 



K„ 



< 00, 



(5.3) 



and there exists an orthonormal basis (ipk)keN of £* such that for all k £ N and all (u,9 n ) G 
H x K n except on a set of potential zer^ 



( Vk ,G n (u:,&^ k ) £ ds (t/ ",e£) = (u,n < lk c{t), 



(5.4) 



where the constants 7fe > 0, independent of n, t and (u,9 n ), satisfy X^eNTfc < °°> an d 
constant C(t) > 0, independent of n, k and (u,0 n ), satisfies lim t _i.oC(^) = 0. 

(D2) For all (3 > and every it holds that 



lim W' 

71— ¥00 



A n (tf;,e?) 



7^|<4., Z "(£)-z"(e ? )) £ |>/3 



M n ((E/T,e?),d£) ds 



o. 



(5.5) 



A set of potential zero is a subset of the state space of the process which the process almost surely never reaches. 
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(D3) Further, for all $ G £ * and all t > it holds that 



lim / E n \($,G(u(s),p(s))$) -a„($,G n (C7 s ",e:)$> | ds = 0. (5.6) 

Finally, we assume that the jump heights of the rescaled martingales are almost surely uni- 
formly bounded, i.e., there exists a constant C < oo such that it holds almost surely for all 
n e N that 

sup v /5^||z n (0?) - z n {@ n t _)\\ £ <C. (5.7) 
t>o 

Then it follows that the process {^Ja n M")t>o converges weakly to an £ -valued centred diffusion 
process characterised by the covariance operator ([5~ 



We now state a second version of the martingale central limit theorem wherein the limiting process 
is characterised by the convergence of the characteristic functions. 

Theorem 5.2. Assume that the laws of the martingales (^a„ M™) t >o form a tight sequence, e.g., 
condition (Dl) is satisfied. 

(D3') The convergence (|5.6[) holds and there exists a sequence (3 n > decreasing to zero such 
that for all $ £ £* 



lim a n E n 

n— ^oo 



A" 

^|<*, Z «(5)-2™(0j)) £ |>/3„ 



<<M n (O-* n (e?)>£|V((£Ce?),d0 d s 



= o. 



Then it follows that the process (y/a n M™) t >o converges weakly to an £ -valued centred diffusion 
process characterised by the covariance operator ([57 



The central condition of the convergence of the quadratic variation processes (|5.6p is unchanged, 
however, the second, technical condition (|5.7[) in (D3) is changed due to the different method of 
proof. That is, condition (|5.8[) arises instead of (|5.7[) as an assumption on the distribution of the 
jump heights employing a characterisation of the limit process using convergence of characteristic 
functions instead of the local martingale problem. The significance for applications of condition 
(|5.8[) in contrast to (|5.7[) is that the former avoids the almost sure uniform bound on the jump 
heights in the latter. That is, arbitrarily large jumps are possible for each martingale in the 
sequence as long as their probability decreases sufficiently fast. Note that (|5.8[) is stronger than 
the similar condition (D2) in the preceding theorem. We omit the proof of the theorem which 
is an adaptation of the estimating procedures in [30J 133 to the infinite-dimensional setting. For 
details we refer to the PhD thesis of one of the present authors [ID] . 

Remark 5.1. We remark without proof that the assumptions (Dl) and (|5.6[) imply the conver- 
gence of the trace processes, i.e., for all T > 

lim a n [ E"Tr G"(t/ S ™, 9") ds = / Tr G(u(s),p{s)) ds . 



5.2 Langevin approximation 

Usually, e.g., in models of excitable membranes, one is ultimately interested in the dynamics 
of the continuous component. We have discussed in Section 12.21 that the coordinate functions 
z™, i — 1, . . . , m, carry all the information needed for the dynamics of the continuous component 
(C/")i>o- Therefore, the knowledge of the coordinate process (z n (Q™)) t >o, or a close approximation 
thereof, is sufficient for many applications. From this point of view the significance of the law of 
large numbers and the martingale central limit theorem is that they provide a justification of an 
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approximation of the processes (J7™, z™(0™)) t >o for large enough n by a deterministic evolution 
equation, on the one hand, and, as we argue in this section, by a stochastic partial differential 
equation on the other hand. 

To this end we first discuss representations of the limiting diffusion in Theorems 15.21 and 15.11 
as a stochastic integral. By definition G(u(s),p(s)) o t is a non-negative, self-adjoint trace class 
operator acting on £, hence there exists a unique non-negative square root, i.e., a non-negative 
operator y/G(u(s),p(s)) o l such that G(u(s),p(s)) 01= y/G(u(s),p(s)) o l o G(u(s) 7 p(s)) o t. 
Let (Wt)t>o be a standard cylindrical Wiener process on £ with covariance operator given by the 
identity (cf. [UGH]). Then, as 

E Tr(y/G(u(a),p( 8 )) ° i\fl) (^/G{u{s),p{s)) o iVl)* ds = / TrG(u(s),p(s )) ds < 00 , 
Jo Jo 

the mapping t > y/G(u(s),p(s)) o i is a valid integrand process for a stochastic integral with 
respect to {W t )t>o- That is, the process (Z t )t>o defined for alH > by 

Z t := / ^G(u(s),p{s)) oLdW s (5.9) 
Jo 

is an £-valucd Gaussian process with continuous sample paths and independent increments which, 
in addition, is also a squarc-intcgrablc martingale. Moreover, the process has the covariance given 
by the operator J G(u(s),p(s)) ds. Therefore, due to unique definition of Gaussian processes via 
their covariance operators, the process (Z t )t>o is a version of the limiting diffusion identified for 
the sequence of martingales ( v / o^'M™) 4 >o. 

Hence, formally inserting the limits into the decomposition of the PDMP we obtain that the 
Langevin approximation (U™, P")t>o of (U", z"(6")) t >o is given by the solution of the system of 
stochastic partial differential equations 

dUp = (A(P t n )U? + B(P t n ,Up))dt 

, (5-10) 
dP? = F(P?,U?)dt+^)jG(U?,P?)dW t . 

The sequence of Langevin approximations (U™,P[ t )t>o possesses the same asymptotic behaviour 
as the sequence of processes (£/", z n (6 n ))t>o- It is obvious that for n —> 00 and thus a n — > 00 the 
noise term in (|5.10[) vanishes and the system approximates the deterministic solution (u(t),p(t)) t >o 
of the system (|4.ip . just as was proven in the law of large numbers Theorem 14 . 1 1 for the sequence of 
PDMPs. It poses no difficulties to make this statement precise in the form of a weak law of large 
numbers similar to Theorem 14. II Thus for large enough a„ one might expect that equation (|5.10l) 
produces a similar behaviour than the PDMP with the major advantage of being analytically (and 
numerically) to a great extent less complex. 

In order to analyse properties of the Langevin approximation, clearly, well-posedness of the sys- 
tem (|5.10p has to be addressed first. This is suitably done within the variational approach to 
stochastic partial differential equations. That is, equation ()5.10[) is assumed to hold as an integral 
equation in X* x £* in contrast to the semigroup approach which defines the solution via the 
semigroup generated by the linear part of (|5.10[) and the variation of constants formula. (Note 
that in its generic form (|5.10[) does not neceassarily posses a fully linear part.) The variational 
approach reflects the approach of using weak solution to abstract evolution equations defining the 
deterministic inter-jump motion of PDMPs taken in this paper. We refer to [3TJ Sec. 1.3.1] for a 
concise introduction to the variational approach to SPDEs containing an existence and unique- 
ness theorem as well as further references. We do not pursue the issue of well-posedness of the 
Langevin approximation any further at this point, as we are of the opinion that this question is 
best addressed when analysing the Langevin approximation for particular models. 

Remark 5.2. The process (|5.9p is not necessarily the only stochastic integral process which is 
a version of the limiting diffusion. Let U be another separable, real Hilbert space, where U = £ 
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is possible, and assume there exists an operator Q £ L\(U, U) (or Q cylindrical) and a function^ 
g £ L 2 ((Q,T),L 2 (U,£)) for all T > such that G(u(i),p(t)) o i = g(u(s),p(s)) o Q o g* (u(s),p(s)) 
for all t>0. Then, the process (zf) t >o defined by the stochastic integral 

zf := / g(u(s),p(s))dW? , (5.11) 
Jo 

where (W t )t>o is an ^-valued Q-Wicncr process, has the same quadratic variation as (Z t )t>o 
hence the processes coincide in distribution. Then starting from the representation (|5.1ip the 
Langevin approximation is given by (|5.10p with the obvious changes in the diffusion term. We 
note that in finite dimensions the non-uniqueness, see, e.g., [3J Chap. 8], of a stochastic inte- 
gral associated to a given covariance matrix can be exploited to improve the speed of numerical 
approximations in Monte-Carlo simulations of diffusion approximations by choosing an optimal 
diffusion coefficient structure, see |32j . In infinite-dimensions the question of a practical implica- 
tion of choosing a diffusion approximation based on (|5.11|) over f|5 . 10[) needs, to the best of our 
knowledge, still to be addressed. 



6 Proofs of the main results 

6.1 Proof of Theorem 14.11 (Law of large numbers) 

The central argument of the subsequent proof is an appropriate application of Gronwall's Lemma 
such that the upper bound satisfies the convergence in probability. Here the estimating procedure 
yielding the estimate to which Gronwall's Lemma is applied necessitates careful attention due 
to more intricate regularity aspects of solutions to abstract evolution equations in contrast to 
solutions of ODEs in Euclidean space. 

The continuous component U" of each PDMP is in between successive jump times the weak 
solution of an abstract evolution equation. Similarly u{i) is the weak solution of the abstract 
evolution equation (|4.1j) . Therefore also the difference of the two paths is in between jump times 
the weak solution of an abstract evolution equation. It thus holds due to jTTJ Sec. 5.9, Thm. 3] 
for almost all t that 

l\\U?-u(t)\\ 2 H = 

2(A(z n (e?)) U? + B(U», z"(©D) ~ A(p(t)) u{t) - B(u(t),p(t)) , U? - u(t)) x . 
Integrating this equation we obtain the integral equation 

K-u(ti)\\ 2 H = K-^Wh (6-1) 
+ 2 r<A(^( ")) U?+B(U?, z n (e:))-A(p(s)) u(s)-B(u(s),p(s)) , U?-u(s)) x ds, 

Jt 

which is valid for almost all to, t% in between two successive jump times. Since both sides of 
equation (|6.1|) are continuous the equality (|6.1I) holds for all to,t\ between successive jump times. 
Moreover, as ?7™ is continuous also at jump times it follows that equation (16. ip holds for all 
t £ [0, T], i.e., we have 

\\U?-u(t)\\ 2 H = \\UZ-u \\ 2 H (6.2) 

+ 2 f(A(z n (e:)) u:+b(u:\ *" (&::)) -a« S )) u( S )-b(u(s),p( S )) , u?-u( 8 )) x d S . 

Jo 

7 Here, L 2 ((0, T), L 2 (U, £)) denotes the space of square-integrable functions on (0, T) taking values in the Hilbcrt- 
space of Hilbert-Schmidt operators from U to £. 



15 



Next we employ the one-sided Lipschitz condition (|4.3[) to estimate the integral in the right hand 
side of equation (|6.2[) . This yields the inequality 



pt ™ pt 

/ \\U?-u(s)\\ 2 H ds + 2L 1 J2 / 
Jo „-_i Jo 



\\UT~Mh < ||Eff-«o|& + 2£i / \\U^-u(s)\\ H ds + 2L 1 ^ / \\z^)- Pj (s)\\ E ds. (6.3) 

3=1 ' 

The overall aim is to apply Gronwall's inequality to the growth inequality (|6.3[) . Therefore, in 
the next step we derive a control on the terms ||z"(9") — Pj(s)|||; in the right hand side of 
inequality (|6.3j) . As p is a solution of (|4.1[) satisfying (|4.2j) we obtain for every functional <f) e E* 
a decomposition 

<<M?(en = (0,^(es)-^-(o)) B (6.4) 

+ [ [A n {<l>,z?(-))E](U?,Q?)d8- f\ < p,F ] (p(s),u( S ))) E ds+(<f>,M?(t)) E , 
Jo Jo 

where the term (<j>, M™(t))E has precisely the form (|3.3|) for all £ 6 [0,T]. Next we expand the 
decomposition (|6.4[) to obtain 



,*?(e?) - Pj (t)) E = (^z?(05) - Pi (o)) B 

+ (<j),Mf(t))E + f ' [A n (<l>,z?(-)) E ](U?,e?) - (^F^i&^^l^Eds 



<0, i^©"), t/™) - *}(p(«), «(»))>* ds . 

We take the supremum over all (j) <E E* with ||</>||.e* < 1 on both sides of this equation, square 
both sides and apply to the right hand side the inequality \a% + ... + a,k\ 2 < fc(|ai| 2 + . . . + |afe| 2 ) 
and the Cauchy-Schwarz inequality which yields 

\\z?m- Pj (t)\\ E 



< 4 ||^(e^)- ft (0)||| + 4 \\M?(t)\\% + 4(^11^(^(6?), U?-FM*)M*))\\ E d*)' 

+ 4(£ || [A n ( ■ , *?(•)>*] if!?, 0") - F 3 (z n (Q^, U?)\\ E ds)" . 
We next apply the Lipschitz condition (|4.4[) on F and obtain the estimate 

\\zj(e?)- Pj (t)\\ E 

pt ™ „t 

< 4 1^(05) - Pj (0)\\ 2 E + 4L 2 / ||C/;- M (.s)|| 2 f d s + 4L 2 ^ / ||<(6?) - K ( S ) III ds 

|[^(.,4X0) B ](^\e; i )-F,(z"(e?),[/: i )|| B d s )\||A/;(t)!||. (6.5) 



To further estimate this term we employ the convergence (|4.5[) of the term ||M™||£; and the 
convergence (|4.6[) of the generator. It follows by the definition of these limits that for every ei > 
and every S > we can find an N ei $ such that for all n > N eiy s it holds due to (|4.5[) for all 
j = 1, . . . , m and all t 6 [0, T] that 

iim;wiu < ^ 

and due to (|4.6[) and the Continuous mapping Theorem that 

(/ T ||[^-,^(0)B](tC,e?)-^-(^(e?),ic)Ld«) a < g 

v./n ' vTt 
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on a set Sli C fl satisfying P n (f2\Sli) < S for all n > N ei) s- Thus continuing to estimate only for 
paths on the set f2i we obtain from (|6.5[) the inequality 



\\zne?)- Pj (t)\\l < 4||z"(e^)-^(o)|| 



5^ 

m 



(6.6) 



4L 



/ \\U?-u(s)\\ 2 H ds + 4L 2 J2 
Jo i=1 Jo 



\\z?W)-Pi(a)\\%d8. 



In order to finally obtain the growth estimate suitable for an application of Gronwall's inequality 
we add inequality (|6.3[) and inequalities (|6.6|) for all j = 1 , . . . , m which yields 



\u? u(t)\\% + J2 \\*?&t)-pj(t)\\% < WK- M 2 h + 4^ ||*?(es) - Pj (o)| 



(6.7) 



3=1 



3=1 



5e 1 + C \\U:-u(s)\\ 2 H ds + C^ 



•Pj(s)|||di 



with constant C = 2Li + AL2m. An application of Gronwall's inequality to (|6.7p yields 



sup - u{t)f H + ]T ||^(e?) - Pj {t)\\ 



te[o,T] 



: ifie 



CT 



3=1 



where 

m 

A"i = |K - u f H + 4^ K( o) - ft(0)lll + 5 £l . 
i=i 

Finally, due to (C3), i.e., the convergence in probability of the initial conditions, it holds that for 
every e-2 > we can find to every 6 > an N £2t s such that on a set f^2 C fi with P n (f2\Q2) < S it 
holds for all n > N e2! s that 

e 2 I, _n,r^n\ ^Ml2 ^ £ 2 



m + 2' 



l (©S)-Pi(o)Hi< 



4(m + 2) 



Vj' = 1, ...,m. 



(6.9) 



Let e, 5 > be arbitrary. Then we obtain choosing e 2 — ee and e\ 
that for all n > 2V e $ := iV ei ,« V iV e2! 5 it holds that 



5(m+2) 



, thus K\ = 62, 



m 

sup te[0 ,T](ll^-^)llff+El!^( e r)-^WIIl) < e 



3=1 

on the set Hi n f^2- Therefore it holds for all n > N e s that 



3=1 

As 5 and e are arbitrary the statement (|4.7[) follows. 



> 6 



25. 



6.2 Proof of Theorem 15.11 (Central limit theorem) 



The proof of Theorem 15.11 is split into three successive steps. In the first step we prove tightness 
of the sequence of martingales which guarantees the existence of a limit. Secondly, we show that 
any limit is a continuous process. Finally, in the last step we prove that the limit is the specific 
diffusion process as stated in the theorem. The conditions (D1)-(D3) in Theorem 15.11 are such 
that each, in addition, to the preceding is needed in the successive steps of the proof. 

Tightness 

In order to prove tightness of the sequence of £ -valued martingales (y / 7i^M™) t >o it suffices to 
show that the following conditions are satisfied, cf. [31] wherein general conditions for tightness of 
sequences of Hilbert space valued processes and, in particular, martingales are considered: 
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(Tl) The sequence of initial conditions (^/ct^M^^o is tight. 
(T2) For all t > it holds that 

lim sup P™ [Tr < v ^"M"> t > S] = . 

5-s-oo ngN 



(6.10) 



and there exists an orthonormal basis (fk)keN °f E* such that for each e > 

0. (6.11) 



lim sup P" [V (ipk^y/a^M^t ip k ) E > e 

m— >oo L^ -- 'k>rn 



(A) The sequence of the real- valued trace processes (Tr <y^M n >()t>o, n > N, satisfies the 
Aldous condition: For every T, e, £ > there exists a h > and an N > such that for any 
sequence of stopping timeCl (c")„>o with a n < T it is valid that 



sup sup P"[|Tr < 1 /a^M n > 

n>N 0<s<h 



-Tr«V^ M ">^ I >S] < e. 



(6.12) 



We next establish the above conditions. First note that condition (Tl) is trivially satisfied as 
Mq = for all n > 0. Hence we proceed to condition (T2). In order to establish the first 
condition (|6.10[) we use Markov's inequality to obtain the estimate 



in [Tr<y^Af Il > t > S] < 



■ E" 



TrG n (r;\ flTl 



(1.S 



where the right hand side is finite due to assumption (|5.3[) . Taking the supremum on both sides 
the same assumption implies (|6.10[) . 

Next, in order to show the second condition (|6.11[) we employ Markov's inequality, the monotone 
convergence theorem (in order to change the order of expectation and the countable summation 
over all k > m), the form of the quadratic variation (|3.11[) and inequality (|5.4[) to obtain for the 
term in the left hand side the estimates 

P"[V { v > k ,<y/^M n » t <p k )s>6] < ^E"[V {<p k ,<M n ^t <Pk)e] 

L ' * k~>m, J n L ' J k*>7n J 



where the upper bound is independent of n G N. Moreover, the property X^gn Ik < oo implies 
that linim^oo J2k>m 7fc = an( i hence (|6.1ip holds for all t > 0. 

Finally, it remains to show (A) . Let T, S > and a n < T be an arbitrary sequence of stopping 
times, then for all for all h > it holds that for s < h 



a„Tr <M n > CTn+s - a„Tr <A/ n > CT „ | > 5 

r<j n +h 



< 



^E"[ 



TvG n (Y r n ,6?)dr 



E 



(<p k ,G n (Y r n ,6?)<p k )edr 



< 



C(h) 



E 

fceN 



Ik- 



Here wc have used Markov's inequality, the strong Markov property of the PDMP and the as- 
sumption (|5.4[) . As the final upper bound is independent of s and n and converges to zero for 
h condition (A) follows. 



'Here every a n is a stopping time on the respective probability space (Q n , T n , (J")t>o,P")- 
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Any Limit is a continuous process 

In the preceding part of the proof we have established that the laws of the sequence of martingales 
(y/ct^, M")t>o are tight which is equivalent to there existence of a weakly convergent subsequence. 
We now prove that under the additional condition (D2) any cluster point of the sequence is 
a measure supported on C(R+,£). The method of proof follows the outline of [23J Lemma 3.2] 
adapted for the stochastic processes being PDMPs on Hilbert spaces, the general setup in this study 
and the particular conditions (Dl) and (D2) in Theorem 15.11 which differ from |23j . Furthermore, 
we have extended the result in |23l Lemma 3.2], which only considers convergence on finite time 
intervals [0, T], to convergence on D(R + ,£). In the following we employ the abbreviations Z™ :— 
yfa^M™ and A t Z n := — Z™_, i.e., (A t Z n ) t >o denotes the process of jump heights. Note that 
A t Z n = ^A t z n (9 n ). 

Further, let P* denote an accumulation point of the sequence (P")„gN- Without loss of generality 
we use P™, n > 1, to also denote the subsequence converging weakly to P*. Furthermore, here 
P™ is understood as a law on the Skorokhod space D(R + ,£) given by the pushforward measure 
of the process (y/a^ M") t > . Then due to the Skorokhod Representation Theorem, e.g., [16j 
Chap. 3,Thm. 1.8], there exists a probability space (f2°, J ro ,P°) supporting D (K + , £) -valued ran- 
dom variables £™, n > 1, and with distributions P" and P* , respectively, such that ( n converges 
to C* almost surely with respect to P°. Further, it clearly holds that E n f(Z n ) = E°/(C") for 
suitable functionals /. 

We begin the proof with preliminary estimates on functions evaluated along the path of the 
PDMPs. These ultimately allow to infer that the process of jumps vanishes in the limit. Let g be 
a measurable, bounded, non-negative function g : M. — > R, that vanishes in a neighbourhood of 
and of oo, that is, there exists a finite constant C' g := swp xeR g(x)/x 2 < oo. For such a function g 
and any $ £ f * we define the process 

GU(*,Z n )s) ■= £ g{($,A s Z n ) £ ) 

se(o,t] 

A"(tc,e?) f ff (v^(*^ n (e)-^(e?)) f )/i ra ((t/r,e™))d S 

' / g(V^($,z n (0 - z n (Q^_)) £ ) M"(d£,ds) , 

JK„ 

where M n is the martingale measure associated with the PDMP, and we infer that G"((<I > , Z n )s) 
is a martingale. Note that the above summation over all s £ (0,t] is well-defined as the PDMPs 
are regular and thus g{(&, A s Z")j) is non-zero for only finitely many s < t. 

The proof now proceeds as follows. We first show (a) that for all t > the random vari- 
ables G™(($, ( n )e), n £ N, are uniformly integrable for all t and (b) that they converge to 
Sse(o t] d{(^> ^sC™)^) m probability. This allows to infer that the convergence result also holds 
as convergence in mean. In part (c) we then use these results to show that the jump heights 
of the canonical process of the law P* arc constantly zero almost surely. This implies that 
F*(C([0, *],£)) = 1 for every t > where C([0,t],£) is understood as the subset of D(R+,£) 
consisting of those cadlag functions which are continuous up to and including time t. The proof 
is completed by (d) extending this result to P*(C(R+,£)) = 1. 

(a) To show that the sequence of random variables G™(($, ( n )e)i n <= N, is uniformly integrable 
in the space (f2°, .F ,P ) it is sufficient that the second moments are uniformly bounded, cf. [THl 
Appendix, Prop. 2.2]. The Ito-isometry for real- valued stochastic integrals with respect to the 
associated martingale measures, which is implied by taking the expectation of the processes in 
22 Prop. 4.5.3], yields 

supE°|Gr(($,C">£)| 2 - S npE n \GU(^,Z n ) £ )\ 2 
= su P E"[/ A"([/?,e?) / g(yfe{$,z n (O-z n (e?))e) 2 vr((U?,0?),a$)d8 
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Therefore, employing the special structure of the map g we obtain the estimate 
supE°|Gr(($,C">£)| 2 

nGN 

< C g sup a r M n \f A n (K,e«) f \($,z n (0-z n (Q^)e\ 2 H n ((U2,Q^,dt)ds 

where the right hand side is finite for every t > due to condition (|5.3[) in (Dl). 

(b) In this part of the proof we establish convergence in probability of the random variables 
G™(($, A( n ) £ ). Let P > be such that g(x) = for \x\ < /?, i.e., the interval (-/3,/3) is contained 
in the neighbourhood of whereon g vanishes. Then we obtain using Markov's inequality and due 
to the boundedness of g the estimates 







h n m,Q n s ) / g ( y fc;(*,z n (z)-z n m) E )n n ((u?,Q«))d8 



< sup a . gR |g(a;)| ^ 



■a^|<#,z«($)- z ™(ey)> E |> l a 



A"([/™,e^) 

Thus due to condition (|5.5|) in (D2) it holds that 

5 (($,A s C ri ) £ )-G t "(($,C") £ ) > <5 



/i"((c/ s ",e«),dc) d s 



lim P°[V 



e(o,t] ' 



= 0. 



Moreover, it holds on (fi°, .F ,P ) almost surely that 

lim ^ g(($,&„C)e) = £ <?«<!>, A S C*>£). 

n— f oo z — * * — * 

se(o,t] «e(o,t] 
Therefore, combining these two convergence results we obtain that 

G»«$,0£)— ► E 3((^A S C> £ ) 
se(o,t] 



(6.13) 



holds as convergence in probability. 

(c) From parts (a) and (b) we infer that (|6.13[) also holds as convergence in mean. Together with 
Jensen's inequality this implies 

lim 



E°(Gr(($,A s C") £ )- S« $ ' A *C*>, 

se(o,T] 

< lim E°|G?«*,A a C n >*)- E 5((^A s C) f ) 

sG(0,t] 



0, 



and hence we infer that 



E° 9({$,KC)e) = lim E°G t "(($,C n ) £ ). 
se(o,t] 



(6.14) 



Furthermore, G™(($, Z n )£) is a martingale which satisfies Gq (($, Z n )s) — 0. This, in turn, implies 
that E n G™(($, = for every n € N. Therefore we obtain due to (|ST4]) 



E* £ g((*,A,Z) e ) E° E ff«*,A a C*>£) 

se(0,t] sG(0,t] 

= lim E°G^(($,C™)£) = lim E n G?(($, Z n ) £ ) = 0. 



(6.15) 
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In a next step, let g m be a sequence of functions satisfying the properties for functions g proposed 
above. Further we assume that the functions g m (x) increase pointwise to x 2 for m — > oo (for an 
example of such functions we refer to |23j ) . Then due to the monotone convergence theorem it 
holds that 

lim E* V g m ((<f>,A s Z) £ ) = E* V \(1>, A S Z) £ \ 2 . 

Tfl — ^ OO 

se(o,t] se(o,t] 

Furthermore, the limiting expectation in the right hand side is zero as each element of the sequence 
of expectations in the left hand side is zero due to (|6.15[) . Next we choose $ to be an element of an 
orthonormal basis (tpk)k<£N of £ and sum the expectations over all elements of the basis yielding 

E E *E se(0 , t] K^z),i 2 . 

fceN 

Due to the dominated convergence theorem we can interchange the countable summation and the 
expectation and, as the PDMP is regular, we afterwards interchange the resulting two summation 
inside the expectation. Then Parseval's identity yields 

e* 1I A ^!H = °- 

sS(0,t] 

As the non-negative random variable inside the expectation is zero only for continuous paths of 
the process (Zs) s e[o,tl we infer that almost all paths are continuous, i.e., P*(C([0, £],£)) = 1. 

(d) To conclude the proof let tk, k € N, be a sequence of times increasing to infinity then 

C(R+,£) = p| C([0,t k ],£), 

and the events in the right hand side satisfy C([0, tfc+i], £ ) Q C([0,tk},£)- The properties of a 
probability measure thus yield 

P*(C(R+,£)) = lim ¥*(C([0,t k ],£)) = 1, 

k— >oo 

that is a process with distribution given by the limit P* possesses almost surely continuous paths. 
Limit is a diffusion process 

In the final part of the proof we uniquely characterise the limit of the sequence of martingales 
{\f(Xn M")t>o under the additional assumptions (D3). The method of proof is via the local mar- 
tingale problem motivated by a proof presented in [34], i.e., the limiting probability measure is 
the unique solution to a particular martingale problem. The author in [33] considers Hilbert space 
valued stochastic integral equations driven by Hilbert space valued martingales with state depen- 
dent quadratic variation. A central limit theorem for the martingales is presented. The arguments 
of the subsequent proof are closely related to [33]. This is as the general result on martingales 
associated with PDMPs, which we have proven in Section [3j result in the problem in this part 
of the proof to be of the same underlying structure as in [33]- One difference, however, is that 
the present conditions (D1)-(D3) are more general than the conditions in [33] and adapted to the 
PDMP setup, hence some estimates differ. 

As in the preceding part of the proof we interpret the sequence of martingales (y/a n M™) t >o de- 
fined on the probability spaces (O n , J-" n , (_F™) t >o, P n ) as random variables on the space D(M.+ ,£) 
equipped with its natural cr-field T>. Further, laws on the canonical space are given by the push- 
forward measure. In order to simplify the notation we denote the laws on the canonical space 
also by P™. Due to results in the preceding two parts of the proof we know the sequence P™, 
neN, admits a limit P* supported on C {&.+ ,£). We use (Ct)t>o to denote the canonical process 
on Z?(R + ,£) which is a version of the martingale (^a n M™) t >o under the push-forward maesure 
P™ for all n £ N or of the weak limit under the measure P* . 



21 



In the following we prove that the limit P* is a solution to a local martingale problem the unique so- 
lution of which is an £ -valued centered diffusion process with covariance operator C(t) G L±(£* , £) 
as given in (|5.2[) . For any twice continuously diffcrcntiable function f : £ —¥ R the extended gen- 
erator Af of such a diffusion is given by 

Af(x,t) = ±Tr(D 2 f(x)aG(t)). 

Then, in order to uniquely characterise the solution to the local martingale problem connected 
with this generator and supported on the space C(M.,£) it suffices to consider mappings / of the 
form (<f>, -)g and ($, -)| for all $ G £*, cf. [31]. That is, we have to show that the canonical process 
ft is such that for all $ G 5* the processes (<E>, and 



<<f,C 



($,G(u s ,p s )$) £ ds 



(6.16) 



are P*-local martingales. We start introducing some notation and then show in parts (a) and (b) 
the local martingale properties of the two indicated processes on the canonical space D(R + ,£). 

As before we use Z" := ^fa^M™ and A t Z n := Z" — Z£_. Further, as indicated above the 
notation is such that we use P™ and E" to denote probabilities and expectations on the original 
given measurable spaces (O n , J 7 ™) as well as on the canonical space (D(M.+ ,£),T>). That is, e.g., 
E n /(Z™) = E"/(£t) for any bounded function /, where the former is the expectation taken on 
the original space (p, n ,J rn ,¥ n ) and the latter the expectation on the canonical space of cadlag 
processes with respect to the pushforward measure. Furthermore, we employ the Ito-formula [33] 
Thm. 25.7] for smooth functions / G C^°(R) applied to semi-martingales. For the particular choice 
of the semi-martingales being the real martingales (4>, Z")s the Ito-formula reads 

f(($,Z?) e ) = 1 y /"(<*, Z s n _) £ )(($,Q„<M n » t *) £ )d S 



s<t 



i^[($,A s z»)^/"(($,zr) 



-M 



(6.17) 



where (M/'™) t > is some martingale on (£2", J 7 ™, (J r ™) t > ,P™) depending on Z n and /. 
Next, we introduce on the canonical space for all positive p the stopping times t p := mi{t G R+ | ||£ t | 
and note that due to the bound (|5.7[) in (D3) on the jump heights we have that for any law 
P™, n > 1, it holds almost surely 

\Kt p \\s<P + C. (6.18) 



Analogously we define the stopping times r" := inf {t G 



\Z?\\ e > p} on the spaces (0", J 7 ™, (F?) t 



Finally, as already mentioned (T> t )t>o denotes the natural filtration on the canonical space. Then 
for A G T> t we define A n := (Z n )~ 1 F G J 7 ™ its preimage with respect to the random variable Z n . 
We now proceed to show that the two processes ($, Ct)f and (|6.16|l are indeed local martingales 
with respect to the limit measure P* . 

(a) Let <E> G £ * be fixed and we choose for every p a smooth function f p G (R) which satisfies 
f p (x) = x if |z| < ||$|| £ .(p + C) and thus /'(z) = 1 and f"(x) = for |z| < \\$\\e-(p + C). 
Therefore it holds for t < t™, which implies the estimate |($, Z^_}s\ < \\3>\\s*(p + C), that 

/;«*, Z?_) e ) = 



and 



/ p (($, Z«) £ ) - / p «d>, > £ ) - ($, A t Z") £ /'(($, Z s "_) £ ) = 0. 
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It follows that applying the Ito-formula (|6.17[) to the function f p and the martingale Z" Al . n all 

terms besides the martingale M n '' p vanish in the the right hand side. Therefore we obtain for 
t<i > t\ and all A £ D tl that 



E" 



l A (($,Ct 2 Ar p ) £ ~ ($,C tl Ar p >£) 



l A (/,(<*, Ct a Ar>) "/,(<*, CtiAr,)*) 



0. 



(6.19) 



The proof of the first martingale property is concluded as in [34]: The mapping £ — ¥ /p(($, Ct 2 Ar }s) 
is almost surely (with respect to the probability P*) continuous and as P™ converges weakly to P* 
it holds due to ([539]) that 



E" 



l A (($,Ct2Ar p >£ - ($,Ct!Ar p >£ 



= 0. 



Here we have employed a weaker version of the continuous mapping theorem, see, e.g., [BJ Thm. 2.7] . 
We infer from the definition of the conditional expectation that the stopped processes are martin- 
gales. Furthermore, as (t possesses continuous paths almost surely under the measure P* it holds 
that t p diverges to oo almost surely for p — > oo. Hence, we can find a sequence of stopping times 
T Pk , k £ N, such that r Ph — > oo almost surely for k — > oo. Thus, the process ($,Ct) is a local 
martingale with respect to P* . 

(b) For the second class of processes we consider smooth functions g p £ C^°(M.) such that g p {x) = 
x 2 for all \x\ < \\&\\ £ *(p + C). Starting from the definition of the conditional expectation as in 
(|6.19[) we obtain 

r-*2AT p 



E" 



1.4 (($,C f2 Ar p ) 



<5>,G(u( S ),p(s))<I>) £ ds~(<Z>Xt^T p y £ 

*,G(u(«),p(s))$> £ . 
,G(u(s),p(s))$) e ds) 



U (S.CfeArJI-^.CnArJI- 



*2AT p 



E" 



tlAT p 

*2AT p 



a n (^,G n (Y s n ,9^) £ ds 



t-iAT p 



E" 



I A n 



a n (^G n (Y s n ,e^) £ - (<f>,G n (u( s ),p( s m £ ds 



tlAr p 



Here the first expectation in the final right hand side vanishes due to the Ito-formula (|6.17p : We 
apply the Ito-formula for the function g p and the martingales Z^ AT n to the terms ($, i^ Al> )| and 

($, Z t " ATp )|. Then we find - similarly to part (a) - that the summands in the right hand side of the 
Ito-formula vanish. Therefore we are left with only the martingale M n ' 9p and the integral term, 
wherein g p ((4>, Z™_) £ ) = 2 for all t < r". The martingale term vanishes due to the martingale 
property and the remaining integral is cancelled by the integral in the above expectation. Overall 
this shows that the first expectation vanishes. 

Next we take the absolute value on both sides of the above equality and obtain, estimating the 
second expectation and extending the integration interval to [0,T], the inequality 



tlAT p 



(H>,G(u(8),p(a))$)eds 



< 



a n ($,G n (y fl n ,^)$) £ - ($,G(u( S ),p( S ))$) £ 



ds. 
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The convergence of the upper bound to zero for n — > oo follows by assumption (|5.6p . Hence we 
have proven an analogous result to (|6.19[) in part (a). The same line of argument that concluded 
part (a) also concludes part (b). The proof is completed. 

7 Application to models of excitable membranes 

The primary motivation for the present work stems from the study of stochastic version of the 
Hodgkin-Huxlcy model [21j describing action potential generation and propagation in spatially 
extended neurons in a PDMP formulation. This model is analogous in structure to hybrid models 
that arc used for the modelling of Calcium dynamics, cf. [IH [24j [42] , or models of cardiac tissue, 
cf. [T5]. Therefore, we consider as an example of the application of the presented limit theorems 
a general compartmental-type hybrid stochastic model for spatially extended excitable membranes 
introduced in gDJ Sec. 3.2] which subsumises the above mentioned applications. (Another example 
for the application of Theorem l4.1l is the law of large numbers that is presented in [5] for a particular 
one-dimensional hybrid model.) We refrain from discussing the physiological derivations of this 
type of model and the implications and interpretations of the limit theorems in this setting. These 
aspects will be subject to a forthcoming publication. We now fix some notation for the remainder 
of the section. The set D C M. d denotes bounded spatial domain with the physically reasonable 
dimensions d < 3. That is, the set D is a bounded interval when d = 1 and when d £ {2, 3} we 
assume it possesses a C 3 -boundary. Further, for a given dimension d, let s denote the smallest 
integer such that s > d/2. Finally, let m £ N be the fixed number of states ion channels can be 
in. 

7.1 Deterministic limit system 

The deterministic limit is the solution to the membrane equation 

d TCI 

il= 2J aij(x)u XtXj + ^gi(x)pi (E t - u) (7.1) 

i,j=l i=l 

with pj , j = 1 . . . , m given by solutions of the coupled equations 

Pj = F j(P, u ) '■= ^Qiji^Pi- Qij(u)Pj ■ (7.2) 

We choose Dirichlet boundary conditions for the component u, i.e., u(t,x) = for all x £ dD 
and all t £ [0, T], which, however, is of no particular importance for the considerations that follow 
and can be readily changed. Here the coefficient functions and are smooth on D, with gi 
non-negative, and the differential operator is strongly elliptic. Further, the rate functions are 
sufficiently smooth|f| Finally, the initial conditions satisfy uo G Hq(D)CiH s (D) andpi(O) £ H S {D) 
and, in addition, the pointwise bounds u(0, x) G u + ] and p;(0, x) € [0, 1], i = 1, . . . , m, such 
that Y^Li x ) — 1) hold for all x € D. Then, the deterministic system (|7.1[) . (|7.2p is well-posed, 
that is, there exists a unique global solution depending continuously on the initial condition, which 
also satisfies g^) EH1 Sec. 3.3.1]. In particular, the solution (u(i),p(*))te[o,T] is in C([0, T],H S (D)) 
componentwise for every T > and pointwise bounded, i.e., u{t, x) £ u + ] and Pi{t, x) £ [0, 1] 
for all (t, x) £ [0, T] x D and alH = 1, . . . , m. 

7.2 Compartmental-type membrane models 

We briefly summarise the essential features of PDMPs (£/",0™) t >o, n £ N, constituting compartmental- 
type membrane models. 

9 In detail the conditions are 1401 Sec. 3.3.1]: The functions qij arc bounded and bounded away from zero on the 
interval [j7_,M4.]. Further, on this interval they satisfy a Lipschitz and polynomial growth condition and are twice 
continuously diffcrentiablc with bounded derivatives. 
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Firstly, an integral component of the sequence of models is a sequence of compartmentalisation of 
the spatial domain D. Thus, for each n G N let V n be a convex partition of the domain D, i.e., 
V n is a finite collection of mutually disjoint conve^ 10 ! subsets of D, called compartments, such that 
their union equals D. 

The second fundamental aspect is the channel distribution across the compartments yielding the 
stochastic jump dynamics and the coordinate functions z n . We assume that each compartment 
either does not contain channels or a fixed deterministic number. Let p(n) denote the number 
of non-empty compartments of the nth model denoted by -Di,n, . . . , D p ( n ),n and l(k,n) be for 
k < p(n) the total number of channels in the fcth non-empty compartment of the nth model. 
Then the piecewise constant components of the PDMPs are given by mp(n)-dimensional vectors 
0™ = (&i' n (t))i=i m ,k=i,...,p(n) with finite state spaces K n . Each component Q^' n (t) counts the 
number of channels located in the domain D^^ n which are in state i at time t. and it holds that 

m 

£e*'"(i) = i(*,n). 
i=i 

as channels can neither be destroyed nor created. We proceed to define the stochastic jump 
dynamics. As two channel switching do not occur simultaneously, the only jumps in the config- 
uration 6 n G K n with non-zero probability are transitions concerning one single channel. That 
is, events for which in one particular compartment one particular channel changes its state. Let 
qij : R — > K + denote the it-dependent instantaneous rate of one channel switching from state i to 
j. Then given a specific configuration 8 n G K n the rate that one channel in compartment Dk >n 
switches from state i to j is given by 

Q§ n (u) eRf, (7.3) 
where Q™ (it) is a functional of the membrane variable it G L 2 (D) defined as 




That is, Q^ n (u) is the instantaneous rate Qij evaluated at the average value of the membrane 
variable over the compartment Dk^ n . Hence the rate (|7.3p is the number of channels in state i 
in domain Dk, n times the rate of one channel switching from i to j. This definition yields by 
summing over all events the total instantaneous rate 

m p(n) 

A>,r) := £ (7.4) 

t,j'=X fc=l 

Note that for each n the total instantaneous rate is bounded and as expected proportional to the 
total number of channels which implies that the PDMPs are regular. Finally, we define on the set 
K n for i = 1, . . . , m the coordinate functions 

P(«) gk,n 

The coordinate process z n (Qf) is cadlag with each component taking values in L 2 (D). Clearly, 
the coordinate process is zero on those compartments which do not contain channels. Moreover, 
each z™(9™) is for every t > a piecewise constant function on the spatial domain D which takes 
values in [0,1]. 

10 The convexity of the compartments is a technical assumption which allows to employ Poincare's inequality in 
the proof of the limit theorems with a known optimal Poincare constant [Tl 1361 . 
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Thirdly, the family of abstract evolution equations (|2.1[) defining the dynamics of the PDMP's 
continuous component U n are given by the parabolic, linear, inhomogeneous second order partial 
differential equations 

d m 

ii=Y, + (Ei - tt) . (7.6) 

i,j=l i—1 

Consistently with the deterministic limit system we equip equation (|7.6I) with Dirichlct boundary 
conditions. Finally, we define the operators A, B depending on 8 n only via suitable coordinate 
functions, cf. (|2.5[) . by 

d m 

A{z n {6 n ))u:= ^ aij (x)u x .... B(z n (e n ),u):^Y,^)^(n(E t -u). (7.7) 

i,j = l i=l 

To conclude, it is easy to see that the characteristics defined via the individual rates (|7.3p . the 
total jump rate (|7.4p and the evolution equation (|7.7[) define a sequence of L 2 (D) x if„-valucd 
infinite-dimensional PDMPs (t/",0™)t>o- Moreover, the membrane component (£/™)t>o is almost 
everywhere pointwise bounded, i.e., U"(x) <E [Ti-,u+] for almost all x £ D and all t > 0, where 
m_ := mirij Ei < and u+ := max^ E{ > 0, for initial conditions [Tq 1 satisfying these bounds, 
cf. [40j Sec. 3.2], which we always assume. 



7.3 Limit theorems for compartmental-type models 

Applying the limit theorems derived in Sections 0] and [5] to compartmental models we find that 
the conditions therein translate into assumptions on the behaviour of the sequence of partitions 
V n and the number of ion channels in the membrane, see Appendix [B] Thus, we denote by S(n) 
the maximal diameter of the non-empty compartments in the nth model, i.e., 

8+(n) := max diam(_Dfc.„) , 

and by t+(n) and I- the maximal and minimal number of channels in non-empty compartments, 
i.e., 

£+(ri) := max l(k,n), £_(n) := min l(k,n) . 

Then the law of large numbers takes the following form. 
Theorem 7.1. Assume that the sequence of partitions satisfies that 

lim 8+(n) = 0, lim t_(n) = oo, (7.8) 

and that the initial conditions (Uq, z u (Qq)) converge in probability to (wo,Po) in the space L 2 (D) m+1 . 
Then the compartmental-type models converge in probability to the deterministic solution of the 
excitable media system (|7.1[) . (|7.2j) in the sense that it holds for all e > that 

m 

j^p[sup te[0 , T] ii^ n -«wiiL» + X! su Pte[o-T]ikr(er)-p(*)ii^ > e ] = o. (7.9) 

i=i 

Moreover, the convergence also holds in the mean in the space L 2 ((0, T), L 2 (D)) ; i.e., 



lim 



|^-«(*)ll^((0 I T),i^)+i;il^( e ")-P(*)ll^((O 1 T),L»)l = 0. (7.10) 
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Next we present the appropriate quadratic variation process for the martingale central limit the- 
orem. For the definition of the limiting diffusion we consider for u,pi € C(D) the bilinear form 



m . 

3=1 ijij JD 

rn „ 

j=i &j Jd 

rn „ 

-Y]Y] Pj(x)q ji (u(x))ip i (x)(t) j (x)dx 

3=1 #3 J ° 
m « 

Pi{ x )^ij{ u { x ))^i{ x )^A x )^ x - 



(7.11) 



3 = 1 vfci 

Note that the right hand side is finite for all 4>i,ipi € L 2 (D) as pt and q_ij{u) are bounded functions. 
Hence, for every given VP £ L 2 (D) m the mapping $ h-> (G(m,p)\E', i>) is a linear, bounded functional 
on L 2 {D) m and, conversely, for every given $ 6 L 2 (D) m the mapping "J M> (G(u,p)^, $) is a 
linear, bounded functional on L 2 {D) m . 

Proposition 7.1. The operator G(u,p) defined via (|7.1ip is /or s > d/2 a trace class operator 
mapping H S (D) into its dual H~ S (D). Moreover, the operator-valued map t H- G(u(t),p(t)) 
defines a unique centred diffusion process on H~ S (D). 

Proof. As stated in [35] it is sufficient for the statement of the proposition that the operator 
G(u(t),p(t)) is self-adjoint, positive and of trace class. These properties are easily verified and for 
a detailed proof we refer to [JO] • d 

In order to state the conditions on the partitions in the central limit theorem we define v+ (n) and 
V— (n) to be the maximum and minimum Lebesgue measure of non-empty compartments, i.e., 

v + (n) := max |-Dfc, n | , V-(n) := min \Dk,n\ ■ 

k— l,...,p(n) k— l,...,p(n) 

Finally, note that in the following the coordinate functions z n arc considered as maps from K n 
into the space H~ 2s (D). 

Theorem 7.2. Let s be the smallest integer such that s > d/2. If in addition to (|7.8[) and the 

convergence of the initial conditions the sequence of partitions satisfies 

£_(n)is-(n) . . 

hm = 1 , (7.12) 

ri-s-oo t + (n) v+yn) 



then the sequence of H 2 " (D) -valued martingales y^J TT^n) ^tj converges weakly to the (H(D) 2s ) m - 
valued diffusion defined by (|7.11[) . 



Remark 7.1. We note that for all reasonable physical domains D and all initial conditions (uo,po) 
sequences of partitions V n an d initial conditions (Uq, Oq) for the PDMPs can be found satisfying 
the conditions of Theorems 17. II and 17.21 For example, a suitable sequence of partitions is obtained 
by grids of uniform cubes with decreasing edge length covering the domain D and putting channels 
only into these cubes which are fully contained in D. For a more detailed discussion of these aspects 
we refer to the PhD thesis of one of the present authors [ID] . 

8 Conclusions 

As a general theoretical results for PDMPs we have derived a law of large numbers and martingale 
central limit theorem in Sections H] and [5] of this study. The former establishes a connection 
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of stochastic hybrid models to deterministic models given, e.g., by systems of partial differential 
equations. Whereas the latter connects the stochastic fluctuations in the hybrid models to diffusion 
processes. As a prerequisite to these limit theorems we carried out a thorough discussion of Hilbert 
space valued martingales associated to the PDMPs. Furthermore, these limit theorems provide the 
basis for a general Langevin approximation to PDMPs, i.e., certain stochastic partial differential 
equations that are expected to be similar in their dynamics to PDMPs. We have applied these 
results to compartmcntal-typc models of spatially extended excitable membranes. Ultimately this 
yields a system of SPDEs which models the internal noise of a biological excitable membrane based 
on a theoretical derivation from exact stochastic hybrid models. 

Topics for further research are motivated by corresponding results in finite-dimensions [3UJ [35] and 
for spatially inhomogeneous chemical reaction systems converging to reaction diffusion equations, 
cf. |25j . In these studies limit theorems are derived for the fluctuations around the deterministic 
limit identified by the law of large numbers. Using the notation of Section[5]we conjecture that the 
sequence of processes, ( V /5^(C/™ — u(t), z"(G") — p(t)) t>0 , n £ N, converges in distribution to a 
suitable diffusion process. Moreover, we further conjecture that this limit is closely related to the 
asymptotic linearisation of the Langevin approximation around the solution of the deterministic 
limit, cf. [33] wherein this result is proven for finite-dimensional PDMPs. 

Further, on the applications side we believe that the Langevin approximation to spatio-temporal 
PDMP models of excitable membranes poses an important object for further investigation. Its 
derivation was the initial motivation of the study of the limit theorems in the present study and it 
is their main application herein which enables to write down the system of SPDEs that constitute 
a Langevin approximation. This system now demands for further analysis, particularly, first of 
all the question of existence and uniqueness of the Langevin approximation has to be addressed. 
Subsequently, as SPDEs arc analytically more accessible than PDMPs a theoretical analysis of 
qualitative and quantitative properties of the models may be possible. 

Finally, we want to mention that the limit theorems presented also find applications beyond 
excitable membrane models. In current work in progress by one of the present authors the limit 
theorems derived in Sections 2] and [5] are applied to stochastic neural field equations, based on a 
model presented in [TJ] , cf. a preliminary account in [55] • We also plan to investigate the connection 
to similar limits derived for reaction-diffusion models, cf. the series of results on variations of the 
model in [25J US [27J [25] and [7J [5J M US ED ■ An answer to this question would contribute to a 
more complete picture of limit-theorems for spatio-temporal stochastic models. 

Acknowledgements: During the time the presented work was accomplished M. Ricdler was a 
PhD student at Hcriot-Watt University supported by the EPSRC grant EP/E03635X/1. M. Riedler 
further acknowledges support from a joint UK Mathematical Neuroscience Network (MNN) and 
the Cell Signalling Network (SIGNET) travel grant. 

A Proof of Theorem 13.11 (Ito-isometry) 

In this proof wc show that under condition (|3.4[) the processes M™, j = 1, . . . , m, n € N, defined in 
(|3.I[) are square-integrable, cadlag martingales which satisfy the Ito-isometry p.5[) . Throughout 
the proof we fix a j = 1, . . . , m and n € N and the results holds for any such j and n. Therefore, 
speaking of a PDMP in the following always refers to the PDMP (J7™,0™)t>o corresponding to 
the fixed n. Further, for notational simplicity wc omit the indices n and j discriminating processes 
and characteristics of PDMPs, i.e., M™ and z™(0") are denoted simply by M and z(Q). Finally, 
recall that t^. k = 1,2, . . ., denotes the sequence of increasing random jump times of the PDMP 
which arc stopping times satisfying lim/ ; _j. 00 Tk ~ oo almost surely. 

First of all, note that the process M is cadlag by definition. The proof of the remaining open results 
is split into three parts. In the first, part (a), we prove the martingale property for the real process 
((<fi, M(t))E)t>o for every <f> € E* , Then, the first main statement of Theorem 13.11 the square- 
intcgrability of the process M(t), is proved in part (b). Moreover, as squarc-integrability implies 
intcgrability, the Hilbert space martingale property follows. Finally, the second main statement, 
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the ItoTsometry (|3.5[) . is established in part (c). The proof we present in part (b) is motivated 
by the proof of [22l Prop. 4.5.3] which states the corresponding results for real- valued martingales 
associated with PDMPs. In extending to the present setup the method of proof employed therein 
one has to ensure, on the one hand, that the employed results and estimation procedures all 
have corresponding analoga in the infinite-dimensional setting. On the other hand, one has to 
carefully make sure that only the weaker regularity results available in infinite-dimensions arc used. 
Finally, the introduction of random initial conditions, not considered in [22], also necessitates some 
adaptations. 

(a) First note that for all <f> £ E* the real-valued processes (</>, M{t)) e satisfy 

(4>,M{t)) E = (<M©*)>i5 - <<M0o))e (A.l) 

- f A([/ S _,6 S _) f {4>,z(Z)) E -{il>,z(Q a -)) Ef i((U a -,Q B -),dt)d8. 

JO JK 

Equation (|A.1I) is obtained from (|3.3[) due to the regularity of the PDMP as the set of jump times 
in [0, t] is almost surely finite for all t > 0. Therefore the integrands in the right hand sides of (|3.3p 
and (|A.1I) differ only on a set of Lcbcsguc measure zero almost surely. Moreover, the integrand in 
the right hand side of (jA.ll) has the form of the extended generator, cf. Theorem 12.11 applied to 
the map 

(«,£)-> (A.2) 

which is independent of u. It follows that the process (</>, M(t))E is a local martingale if the map 
(|A.2j) is in the domain of the extended generator, cf. Thcorcm l2.ll Obviously, path-diffcrcntiability 
almost everywhere is trivially satisfied as the map t n- ((f), z(@ t ))E is piecewise constant. Hence, 
it remains to consider the integrability condition for which it is a sufficient that 

E I A(U a -,Q s -) I \(<f>,z(£) — z(@ s -))e\ (J-((U s -,Q s -),d£)ds < oc Vt>0, (A.3) 
Jo JK 

cf. [131 [HI [22] . Using Young's inequality we obtain an upper bound to (|A.3[) by 

Ie[ a([/ s _, e s _) ds + W A(t/ S _, e s _) / 1 (0, z(o - Z (e s -)) E \ 2 K(u s -,o s -),d0 ds. 

1 JO 1 JO JK 

Here the first expectation is finite due to the PDMP being regular and the second is finite by an 
immediate consequence of assumption p. 41) . 

Next, wc show that the process is not only a local martingale but even a martingale. As mentioned 
above the process (</>, M (t)) e satisfies 

(<t>,M{t)) E = [ I (&z(fl-«(e._)) JS M(ds > dO 
Jo JK 

where M := N — N is the random martingale measure associated with the PDMP with count- 
ing measure N and compensator Af(d£,ds) = A([/ S _,0 ;j _) , O s -), d£) ds. The validity of 
this formula follows as ()A.2[) is in the extended generator. Thus the process (4>,M(t)} e has the 
form of a stochastic integral with respect to the martingale measure associated with the PDMP. 
Furthermore, due to [2H Thm. 4.6.1] it holds that the process is a martingale if (|A.3j) is finite 
for all i > 0. But we have already shown that this holds due to the regularity of the PDMP and 
assumption (|3.4[) . 

(b) We now prove the square-integrability of the process M. In a first step we prove in (b.l) 
that M stopped at the first jump Ti is square-intcgrable. Subsequently in part (b.2) this result is 
extended to M stopped at any jump time r^, k G N. Then we are able to infer square-integrability 
of the process M. As square-integrability implies integrability it follows from part (a) that M is 
a Hilbert space valued martingale. 
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(b.l) Note that prior to t\ the jump component of the PDMP remains constant. We introduce 
the notation 



A(U r ,Q ) / z(0-z(6 )n((U r ,e ),d0dr 



N(s) 



which implies that s i— > ||-^V(s)lH i s almost surely absolutely continuous with derivative 



Ts \\N(s)\\% 



■ dt 



N(s),N(s)) ] 



2A(U s ,e ) / (z(t)-z(Q ),N(8)) E ii((U a ,e a ),d$). 



(A.4) 



K 



Due to the structure of a PDMP we obtain for the conditional expectation with respect to the 
initial condition 



E 



[||M(riAt)||||Jb] = \\N(t)\\l exp(- t A([/ r ,9 ) 



i dr 



K 



\z^)-z(e Q )-N(s)\\ 2 B ^((U s ,Q ),d'&)) A([/ S ,e )exp(- / A(U r , 6 ) dr ) ds 



That is, the first term in the right hand side is the position of the stopped process ||M(ti A t)\\% 
at time t if t < t\ times the conditional probability that the first jump does not occur before t. 
The second term is its position after the jump integrated over the conditional density that a jump 
occurs in [0,i]. We apply integration by parts to the hrst term (note that N(0) = 0) and find that 



\\N(t)\\% exp - / A(U Si e )ds 



2(£N(s),N(s)) E exp - / A([/ r ,e )dr 



N(s)\\%A(U a ,e ) exp 



A(C/ r ,6 )dr 



ds. 



Therefore we obtain 

E[||M(nAt)||||Jb] 



2(±N(s),N(s)) B exp(- / A(U r ,Q )dr 



ds 



A 



- z(Q ) - N(s)\\l - \\N(s)\\ 2 E t,((U s ,Q ), &&)) 
A(U S ,Q ) expf- / A(C7 r ,9 )dr 



Note that - z(6 ) - N(s)\\ 2 E 

thus 

E[\\M(nAt)\\l\F ] = 



ds. 



|z(0) - z(0 o )||| + - 2[z{&) - z(Q ), N(s)) E and 



/' 

Jo 



2(±N(s),N(s)) E exp[- 



A(C/ r ,6 )dr) 



/ r 



A' 



A" 



>(0)-z(e o ), N(s)) E fi((U s , O ), d0))J A([/ s , 0„) exp^y A(l7 r , 6 ) dr 
|z(0) - «(e )||| A*((^», @o), dt?))) A(f/ S , Go) expf- f A(U r , O ) dr 



ds 



ds. 



Due to form of the derivative (|A.4[) the first two terms cancel and we are left with the equality 
E [|| M( n At) || || F ] (A.5) 



A(U S ,Q Q ) f \\z(#)-z(Q )\\ E u((U s ,e ),dti)) exp(- f A(U r , 6 ) dr) 



Ids. 
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Next we calculate the expectation of the real- valued process 



N 2 (s) 



A(U r ,Q ) / \\z(#)-z(Q )\\%(i((U r ,Q ),dti)dr 



K 



stopped at t\. The process N 2 is connected to the process N defined at the beginning of part 
(b.l) inasmuch as the integrand of the former is the squared norm of the latter. Furthermore 
note that N 2 is the term inside the expectation in the right hand side of the Ito-isometry (|3.5[) . 
Thus the aim is now to show that the conditional expectation of N 2 (t A Ti ) equals the conditional 
expectation of || M{t A Ti)||f;. Again due to the particular structure of the PDMP we obtain for 
the conditional expectation 

E[7V 2 (n At) I To] 

= N 2 (s) cxp(-^ A{U r , e ) dr) +^ [iV 2 (s)A(C/ s ,e ) cx p(-^ A([/ r ,6 )dr) 
Integration by parts applied to the integral term yields 

N 2 (s)A(U s ,Q ) ex p(-^ A(U r ,6 )dr)] ds = -iV 2 (t) exp(-^ A(U r ,Q )dr 
A(U a ,e ) [ \\z(0)-z(e o )\\%li((U B ,Qo),d0) exp(- ^ A([/ r ,e )dr) 



ds. 



d.s. 



Therefore we obtain that 
E[N 2 (nAt)\T ] 



(A.6) 



A(U s ,e ) Jjz^)- z(Q )\\%fi((U s> Q ),d-d) exp(-^ A(U r , 9 ) dr) 



A comparison of the right hand sides in equalities (|A.5[) and (|A.6[) shows that they are equal and 
thus we obtain after taking the expectation of both conditional expectations that 



|M(nAt)|| 



lN 2 (n At) 



(A.7) 



As A^ is increasing and thus A*2 (ti At) < N 2 (t) almost surely, we obtain that the right hand side 
in this equation is finite due to condition (|3.4p . Note that (|A.7j) is the Ito-isometry p. 51) for the 
stopped process M{t A t{). 

(b.2) In this part of the proof we show the square-integrability for the process M stopped at an 
arbitrary jump time Tfc, k £ N, and finally for the non-stopped process M. To this end we first 
note that Analogously to part (b.l) we find that 



E 



\M{T k+1 At.)- M(T k At)\ 



E I ^ T k 



= E 



N 2 {r k+1 At) - N 2 {r k At)\F T 



Thus taking expectations on both sides of this equality yields 

, 2 



E\\M{ Tk+1 At)- M{r k At) 



E 



N 2 (r k+1 A t) — E N 2 (r k A t) < oo 



(A., 



where the right hand side is finite as due to (|3.4p both expectations are finite. 
By induction we next show that each M{r k At) is square-intcgrablc. Assume that E||M(r/ £ Ai)|||. < 
oo, where the induction basis for k = 1 holds due to part (b.l). Then the reverse triangle inequality 
yields that 

E||M(t*+i A t)\\% + E\\M(r k A t)\\% - 2E(||M(r fe+1 A t)\\ E \\M(r k A t)\\ E ) 



< E\\M{T k+1 At) - M(r k At)\\] 
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Here the right hand side is finite due to (|A.8|) and an application of Young's inequality to the 
product in the left hand side yields that for all e > 

(l-2 e )E||M(r fe+1 At)|||+(l-i)E||M(r fc At)||| < 00 

Assume that E||M(rfc+i At)|||; = 00. Then choosing e < 1/2 we obtain a contradiction due to the 
induction hypotheses. 

In a final step of this part of the proof we show square-intcgrability for the non-stopped process. 
Using Fatou's Lemma and monotone convergence for interchanging limits and expectation we 
obtain the following upper estimate 

E||M(f)||| = Eliminf ||M(r fe At)||| 

k— ¥00 

< liminf E\\M(r k At)\\% = lim EN 2 {r k A t) = EN 2 (t) , (A.9) 

k— >oo k— >oo 

where the final term is finite due to condition (|3.4p . Moreover, as square- integrability implies 
integrability, the martingale property for the Hilbert space valued process M now follows due to 
part (a). 

(c) Finally in the last part of the proof we establish the Ito-isometry. To this end we first show 
that equality (|A.7|) holds for all r k A t, k <E N. Again we proceed by induction with the induction 
basis given by (|A.7[) . We observe that 

\\M{T k+1 A t)-M(r k A t) HI = \\M{T k+1 At)\\ 2 E -\\M(T k At)\\% 

- 2(M{r k A t),M(r k+1 At) — M(r k A t)) E . (A.10) 

Taking the conditional expectation with respect to the stopped fi-field F Tk At we nn d that the 
second term in the right hand side of (|A.10[) vanishes as it holds 

E[(M(i* A t), M(r k+1 At)- M(r k At)) E | T TkAt ] = 

= (M{T k At),E[M{T k+1 At)-M(T k At))\F TkM ]) E = 

due to the following properties of the conditional expectation: Firstly, for E-valued random 
variables X,Y such that E||X||£;||Y||.e < 00 it holds for ^/-measurable X that E[(X, F)e|0] = 
(X,E[Y'|5]) £ , [HJ Lemma 2.1.2]. Secondly, the Optional Sampling Theorem, i.e., E[M(r k+ i A 
t) J JvfcAt] = M{ T kAt)) in the above application, also holds for Hilbert space- valued martingales^- 
Thus we obtain 

E[||M(r fc+1 A t)\\ 2 E \T TkM ] -E[\\M(r k A t)||||^ A t] 

= E[N 2 {T k+1 A t)\F TkM ] - E[N 2 (r k A t)\F TkAt ] . 

Taking the expectation on both sides of this equality and using the induction hypotheses, i.e., the 
second expectations on both sides of the above equality equate, yields 

E||M(t*+i A t)\\% = EN 2 {r k+1 A t) . (A.ll) 

We conclude the proof extending the Ito-isometry (|A.11[) from the stopped processes to the non- 
stopped process. We have already obtained the upper estimate E||M(f)|||; < EN 2 (t), cf. (|A.9|) . 
Hence it remains to prove that a lower bound is given by the same term. As ||Af(£)||^ is a real- 
valued submartingale it holds for all k > 1 due to the standard Optional Sampling Theorem for 
cadlag submartingales, see, e.g., [55] App. B], that 

E\\M(t)\\ 2 E > E||M(r fc At)||| = EN 2 (r k At) . 

Hence, for k — > 00 we obtain by monotone convergence E||M(f)|||; > EA^^) which, combined 
with the upper bound (|A.9j) . yields the Ito-isometry (|3.5[l . The proof is completed. 

11 The Optional Sampling Theorem can be proved similarly to the methods employed for 1411 Lemma 2.1.2] relying 
on the linearity properties of the Bochncr integral and the monotone convergence theorem. 
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B Proofs for the neuron models 

B.l Proof of Theorem HJ] (Conditions for the LLN) 

We apply Theorem 14. II for the choice of spaces X = H^(D), H = L 2 (D) and E = L 2 (D). Hence, 
we have to prove in the following that the assumptions therein are satisfied, i.e., (i) the one-sided 
Lipschitz condition (|4.3|) on the operators A and B defined by (|7.7p , (ii) the Lipschitz condition on 
the right hand side of the gating system (|7.2j) , (iii) the uniform convergence of the generator and 
(iv) the martingale convergence. Finally, in (v) we extend the convergence in probability due to 
Theorem 14.11 to convergence in the mean (|7.10[) . In the following we use • to denote the pointwise 
product of real functions on D. 

(i) For the non-linear operator B we find that the left hand side in the Lipschitz condition is for 
almost all t given by a finite sum of terms 

(Pi- (Ei-u) -pi- {Ei-v),u-v) H i , i = l, ...,m, (B.l) 

with u, v G Hq(D) and Pi, Pi £ L 2 (D). Hence, the duality pairing corresponds to the inner 
product in L 2 (D). We estimate each of the summand of the type (IB. II) separately. Using the 
triangle inequality we obtain 

\(pi ■ (Ej - u) - -pi ■ (Ei - v),u - v) H i\ < \Ei\ \ (pi-Pi,u-v) L i\ + \(pi-u-pi- v,u -v) L 2\. 

Here, the first term in this right hand side is further estimated using Cauchy-Schwarz and Young's 
inequality, which yields 

\(pi -pi,u-v) L 2\ < l\\pi-Pi\\ L2 + l\\u-v\\ L2 . 

For the second term we obtain, making use of the triangle inequality, Cauchy-Schwarz and Young's 
inequality and the pointwise bounds on pi and i>, the sequence of estimates 

\(pi-u-pi-v,u-v) L2 \ < \{pi ■ (u - v),u - v) L 2\ + -pi,v ■ (u - v)) L 2\ 

< \\pi ■ (u - v)\\ L2 \\u -v\\ L3 + \\pi -pi\\ L2 \\v ■ (u - v)\\ L2 

II II 2 i7 2 ll II 2 111 ~ll 2 

- \\ U ~ V \\l2 + T \\ U ~ + 2\\Pi ~PA\li ■ 

A summation over all these estimates for i = 1, . . . m yields 

(B(p,u) - B{p,v),u-v) H i < m(l + 2^) ||u-«||i a + \\ Pi ~P*\\l» ■ 

i=i 

Adding the estimate 

(A(u-v),u-v) H i < -Ji\\u - v\\ 2 H i + j 2 \\u - v\\ 2 L 2 < J2\\u - v\\ 2 L 2 

for some 71,72 > 0, which holds as the linear operator A is coercive and independent of p, we 
obtain 

m 

(A(u - v),u- v) H i + (B(p,u) - B(p,v),u- v) H i < C (||u - v\\%2 + ^ \\pi -PtH^a) 

i=l 

for a suitable constant C . Finally, integrating over (0, T) we find the one-sided Lipschitz condition 
(|4.3p is satisfied. 

(ii) Due to the triangle inequality it suffices to consider differences of the form \\pi ■ q(u) — pi ■ 
3(l>)||z,2, where q substitutes for an arbitrary rate function qjk- Using the triangle inequality, the 



1 1 — m 

1 + U sr-^ 11 ^ ||2 
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pointwise boundedness of pi and q by 1 and q, respectively, and the Lipschitz condition on the 
rate functions q (with common Lipschitz constant L) we obtain 

\\Pi ■ q(u) - Pi ■ q(v)\\ L 2 < \\pi ■ q{u) - pi ■ q(u)\\ L 2 + \\pi ■ q(u) - p, ■ q{v)\\ L 2 

< q \\pi -Pi\\L 2 + l \\u - v\\ L 2 . 

A summation over all such separate estimates, integrating and squaring both resulting sides yield 
the Lipschitz condition (|4.4[) . 

(iii) In order to prove the convergence of the generators (|4.6p we employ in the following two tech- 
nical results which we collect in a separate proposition. Firstly, the purpose of the formula (|B.2|) is 
to transform the generator of the PDMP into a form that allows comparison with the deterministic 
limit system (|7.2p . Secondly, the inequality (|B.3[) . which bounds the norm |j?7™|iL 2 ((o,T),ff 1 ) by 
a deterministic constant uniformly over n G N, is used repeatedly in the subsequent estimation 
procedures. 

Proposition B.l. (a) The generator of the PDMP satisfies 




where 

p(n) 

fc=l 

(b) For all n G N and all T > it holds that 

/ \\U t n \\ 2 Hl dt <d(l + T)e 2C2T , (B.3) 
Jo 

where the constants C\^G% are deterministic and independent of n G N. 

Proof, (a) We denote by 9^ i-v ■ for all k = 1, . . . ,p(n) and all i ^ j, i,j = 1, . . . m the config- 
uration in K n that arises from the configuration 9" through the event that a channel in state i 
located in the compartment Dk, n switches to state j. Then simple reorganisation of finite sums 
yields 

p{n) p{n) 

= EE(^"(W - °f n + E E (*rosw - *r(n) # n 

fe=l j^i fc=l j^i 

- EE(^'»...) <• <*•(-) + EEt-jj^ 1 ^) «5» 

fc=l j^i V ' y fc=l V ' ; 

P(") P(") 

= E^cn^E^rw 1 ^,-) -I>r(n-(E^>)i^,„)- 

j^i fe=l j^i fe=l 

Thus we obtain that the generator satisfies (|B.2I) . 

(b) By definition of a PDMP it holds that the component (U")t>o is the weak solution of the 
evolution equation 

m 

u? = au? + ^> *?(e?) - U?) 

i=l 
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with initial condition Uq. We consider the reaction term in this equation as a given inhomogene- 
ity. Then standard estimation procedures from the theory of linear parabolic partial differential 
equations, cf. [171 Sec. 7], yield, after appropriately estimating the inhomogeneous term, 

£ \\u;Tm & < K lC ^ T (\\uz\\h + 2« 2 £ IMU~ £ |K(e?)|Ui dt) , 

where the constants K\ , K2 are deterministic and depend only on the domain D and the coefficients 
of A. Further, it holds that < and the sequence of initial conditions is bounded 

by assumption as Uq(x) € [u_,u^.] for all x £ D almost surely. The inequality (|B.3[) follows. □ 



We now proceed to the actual proof of the convergence (|4.6D . To this end we need to consider for 
almost every t and all i = 1, . . . , m, the convergence in L 2 (D) of (|B.2|) to Fi(z n (Qf), f/™) where 
Fi is as defined in (|7.2p . That is, we have to estimate 



E 



z«(er) • q -(un - *?(e?) • (c/f )) - E ( z "( e ") • fciW) - ^"( ") • WW) 



Jt-> 



(B.4) 



We find that the single summands in the two summations match up and thus it suffices to 
consider each of them separately. Employing the boundedness of the coordinate functions, i.e., 
||z"(8")||_loo < 1 we obtain the estimates 



\zne?)-q%(ur)-zne?)- qji (u? 



\L 2 



< 



p(n) 



mm 



k=l 



Qii{Ut)\\ 

-quun 



L 2 
2 

L- 



p{n) 

E 

k=l 



Q^m-liAU^Wdx. (B.5) 



For the last equality we have used that the summands are mutually orthogonal in L 2 (D). Next 
we estimate each of the remaining integrals in (|B.5|) using the Lipschitz continuity of qtj and 
Poincare's inequality in L 2 (Dk t n), i.e., 



D, 



<1>J 



\Du n \ 



D k . 



V?(v)dv)- qij (U?(x)) 



d.r 



< L 



\D 



k,n\ JD 



U?{y)dy-U?{x) 



d.r 



< i^" 2 diamp fc . n ) 2 ||V[/; l ||i 2( ^ n) , 

where ||VJ7"||ia/£) fc \ is the norm in L 2 (Dk n ) of the Euclidean norm of the gradient vector 
VJ7". Here we have employed that for convex domains the optimal Poincare constant is given 
by 7r _1 diam(Z?fc „) |36j . Hence, a summation over all k = 1, . . . ,p(n) and employing the estimate 

livc/f ||| 2 < \\uhh yields 

\\q?i(U?) - qji(V?)\\ 2 La < 6 + (n) 2 L 2 Tr- 2 \\Un\ 2 m ■ 
Integrating over (0,T) we therefore obtain for (|B.4p the estimate 

rT 

U(d>, z?(-))r.*](u?,e?) - FAz n (e?\ur) 

L' 2 

\2 r2„-2 







dt 



< 5 + (n) 2 L 2 ir- 2 2(m - 1) / || U?\\%i dt . 
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Finally, the norm ||C/™||L 2 ((o,T),_f/ 1 ) is bounded independently of n £ N by a deterministic constant 
due to Proposition lB.lf b'). This upper bound holds for almost all paths of the PDMPs (J7™, 8")t>o 
and thus there exists a constant C > independent of n such that 



[A^z^^iU^e^-F^ie^Un L2 dt < 6 + (nfC 



(B.6) 



almost surely. Due to the assumption (|7.8[) the estimate in the right hand side converges to zero 
for n — ► co and the convergence (|4.6p follows. 



(iv) Next we consider convergence in probability of the martingale part. To this end we employ 
Lemma 13.11 As before we denote by 9^ the channel configuration that arises from the con- 
figuration 8 n if a channel in compartment Dfe in switches from state i to state j. Then it holds 
that 



p(n) 

= E E(l W^;( s )) - z i(Qs)\\i* Q-f(K) e?(s) 
k=i j^i 



p(n) 

* *E 



1 5>?(a) + eyw) 



Z(fc,n) 2 
fe=i v ; j#t 



This implies that 



A*(EC,e?) / \\zm-4m\\h^{(K^),^) &s = 0(l-{n)- x ) 



K„ 



Hence, under condition 



the assumption of Lemma 13.11 is satisfied. 



(v) Finally, we extend the convergence in probability to convergence in the mean for the indi- 
vidual components being in the space L 2 ((0, T), L 2 ), see the remark following Theorem 14. II First 
of all note that the components arc bounded, i.e., 



\\U?-u(t)\\ L z<2u\D\, 
Therefore it holds that 



\\z?(®?)-pM\\l*<2\D\ 



\x n -x\\ ;=\\u n -u\\ L2mTlL 2 )+ J2 lk"(0") -j*IUw),i») ^ c 



for a suitable deterministic bound C < oo independent of n £ N. Then for all eo > it holds that 

E n |LY„-A|| = E^II^-XIII^^H^]] +E"[||A„-A||I [ || Xn _ X || >eo] ] 

< e +MP™[||X„-jq >e ] ■ 

Next choose eo < e /2 and note that due to the convergence in probability there exists an N t such 
that MP™ [\\X„ - X\\ > eo] < e/2 for all n > N e . Hence, for every e > there exists an N € such 
that E"||A„ — A || < e for all n > N e . Convergence in the mean is proven. 

B.2 Proof of Theorem HH (Conditions for the CLT) 

In order to prove Theorem 17.21 we employ Theorem 15.11 for the space E = H~ 2a (D) where s is 
the smallest integer such that s > d/2. We usually employ the simpler notation E and £ = 
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E m throughout the proof, however occasionally switch to H~ 2s (D) if we want to emphasise the 
specific choice of the Hilbert space. The reason choosing this particular integer s is that it is 
the smallest integer such that the embedding of H 2s (D) into H S (D) is of Hilbcrt-Schmidt tvptf^l 
due to Maurin's Theorem [31 Thm. 6.61] and H S (D) is embedded in C(D) due to the Sobolev 
Embedding Theorem. These two properties are essential in order to prove the conditions (|5.3[) 
- (|5.7[) of Theorem 15.11 All conditions except (|5.6D , which establishes the convergence of the 
quadratic variation, arc straightforward consequences of the assumptions of the theorem. These 
are shown in part (i) of the subsequent proof. For condition (|5.6I) more involved estimation 
procedures are necessary which are presented in part (ii). 

(i) We first show condition (|5.3[) . As in the preceding section @ k i->-j denotes the element of K n 
that differs from 6 n by one channel in the kth. compartment being in state i instead of state j. 
Then, the Sobolev Embedding Theorem yields the estimate 

\\*m^M) ~ z n (®?)U = sup iKMr^WH-l <77^—,\D k > n \, (B.7) 



U\\ H 2s 



l(k, n) 



where C is a constant resulting from the continuous embedding of H 2s (D) into C(D). Using this 
estimate for the jump heights in the space H~ 2s (D) we find similarly to part (iv) of the proof of 
Theorem 17. II that it holds 

\k n {ui\&i) j ||z"(0-^(eni/i"((c/; i ,er),dOdil =0(1). 

Hence, condition (|5.3[) is satisfied. Moreover, we infer from (|B.7p that the rescaled jump sizes 
are bounded almost surely uniformly, i.e., condition (|5.7p is satisfied. Particularly, it holds that 
y/a^\\z n {0k,i->j) ~ z n {.9 n )\\ E = Olefin)- 1 / 2 ). This implies that for arbitrary /3 > and any 
$ g (H 2s (D)) m there exists N p such that for all n > N p 

^((u.n.d*) = 

W|(*,« 4 n (£)-*?(« n )>el>0 

holds for all values (u,8 n ) the PDMP attains. Therefore, by dominated convergence we infer 
that also condition (|5.5|) is satisfied. It remains to consider condition (|5.4[) . To this end let 
(Vfe)fcGN be an orthonormal basis in (H 2s (D)) m , where ipk = {fl, ■ ■ ■ , P™) and hence (ip l k )keN is 
an orthonormal basis in H 2s (D) for alH = 1, . . . , m. Then we obtain the estimate 

(<p k ,G n (u?,e?)< Pk )s 

= a"([/« en / (Y,{vizm-^)) H -A v n m\&i i )A0 

m p 

< mY,M.\\H°U n (u?,&?) / \\z?®-zne?)\\ 2 H-.i* n {{u?,W),<%) 
i=l v Jli ~ 

Here we have employed for the the individual summands in the right hand side that for zf(£) — 
z™(9") g L 2 {D) the duality pairing in H 2s {D) equals the duality pairing in H S {D). Further, note 
that (£) — z^(6^)\\u-s satisfies an estimate analogous to (|B.7|) due to the continuous embedding 
of H S {D) in C(D). Therefore we overall obtain that 

m 
i=l 

for a suitable non-random constant C independent of n. Finally, set jk X)i=i ll^fclllf 3 then it 
holds that X^fceN 7fc < oo as the embedding H 2s (D) H S (D) is of Hilbert-Schmidt type. We 
infer that condition (15.41) is satisfied. 



1 The embedding of a Hilbert space X into another Hilbert space H is of Hilbert-Schmidt type if JZfegN 
oo for every orthonormal basis (i/?fc)fcgN of X. 
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(ii) In the second part of the proof we establish the central condition (|5.6[) of the convergence 
of the quadratic variation. For simplicity of notation we omit the time argument of the PDMP 
paths and the deterministic solution as the following estimates hold for almost all t. First of all 
we expand the quadratic variation of the martingales into the finite sum 

A n (u n ,e n ) [ ($,z n (0 -z™(e n ))i// l (([/ n ,e"),d£) 



m m p(n) Qk,n 



m m p(n) p^k.n 

=£ E E mfa oSr^Wi.ioi + EEE Q-r(u n )^j Dk j 

3=1 i=i k=l v ' ' 3=1 i— 1 fc=l v ' ; 



2 



(£/«) 



B 



k.n 



Q k u n (u n ))(^D k jE(4> 3 ,i Dk J 1 



l(k, n) 



2 



We find that the terms in this summation match with the integral terms in the definition of 
the operator G(u,p) in (|7.11[) . Thus, due to the triangle inequality it suffices to consider the 
convergence of the single summands separately, i.e., we have to consider, on the one hand, for all 
j = 1, ... ,m and i ^ j the differences 



D 



p(n) Qk,; 

Pj(x)qji(u(x))4>j(x)dx - an/. 



k=l 



l{k, n] 



2 



l (u n )(4> 3 J Dk J 



and, on the other hand, for all i, j = 1, . . . , m such that i ^ j the differences 



D 



Pi(x) qij{u) </>i(x)<l>j(x) dx - a n V 

k=i L ^ n > 



(Bi 



(B.9) 



We next estimate these terms separately in parts (ii. 1 ) and (ii.2). Finally, in part (ii.3) the 
estimates are combined to prove the convergence of the quadratic variation, 
(ii.l) A further application of the triangle inequality yields 



KB = 



Pj(x) qji(u(x)) (t>]{x) Ax 



3 



+ lgl(ifc) J Dk J^x))^x)dx - Ong^QjiWi.I 



D k , n )E 



(B.10) 



We estimate the two resulting differences separately and obtain for the first term in the right hand 
side of ([B.lOp the estimate 

p j (x)q ji (u(x))<p 2 j (x)dx - f z"(Q n )(x)q Jt (U n (x))0 2 (x)dx 



< 



p j {x)q ji {u{x))(l) 2 j (x)dx- / z™(O n )(x)qji(u(x))(j) 2 (x)dx 



D 



D 



z^{@ n ){x)q ji {u{x))4> 2 {x)dx- / z^{@ n ){x) qji {U n {x)) 4> 2 {x) dx 



D 



< n4>i\\l~\\vi-zU& n )\W + LUi\\l-h-u n \ 



L 1 



(B.ll) 



For the second term in the right hand side of (|B.10[) we obtain by employing Qj' n /l(k, n) < 1 the 
estimate 

P( n ) n. 

E / qji{u n (x))^{x)ax 



k=l 



l{k, n 



38 



and we continue estimating each summand therein separately. We begin employing the Mean 
Value Theorem to expand the rate function qji in the integral in the left hand side such that 

qji (U n (x))= qij (-}— [ U n (y)dy)+q' ji (# k > n (x))(u n (x)--j— [ U n (y)dy), (B.13) 

where i? fc '"(x) denotes an appropriate mean value. For now we omit the remainder term, i.e., the 
second term in the right hand side of (|B.13j) . a consideration of which is deferred. Hence, we 
obtain for the absolute value in each summand in ()B.12|) the estimate 



ir^—, [ U n (y)dy) f <j>%x)dx 



l(k,n) \J Dkn 

We note that qji is bounded by q and continue estimating which yields 

1 



4>j (x) dx 



< q\D k , 

< q 



q\D k 



\D k ,n\ 

\ - 



D, 



1 



an\Dk,n\ 2 

l(k,n)\D k>n \ 



D 



k,n\ JD 



4>j (x) dx 



(B.14) 
(B.15) 



The term (|B.14|) is estimated using Poincare's inequality which yields an upper bound by 

n over all k = 

£-(n) V-{n) 



qw 2 diam 2 (Dfc,n) || V</>|| 2 2 ( Dk y For the terms (|B.15|) a summation over all k = 1, . . . ,p(n) yields 



p(n) 

q E \ D >« 



k=l 



1 



a n \D k J 2 



l(k,n)\D k , n \ 



D 



k,n\ JD 



4>j (x) dx) < q 



1 



l + (n)v + (n) 



(B.16) 



where </>" is a piecewise constant approximation to <pj defined by 

p(n) 



(\D k , n \ J D 




4>j (x) dx^j 







As (pj converges to <f)j in L 2 (D) it holds that the sequence of norms converge, hence ||<^>™||l 2 is a 
bounded sequence. Therefore the right hand side in (|B.16I) is a componentwise product of con- 
vergent sequences. The sequence |1 — (£-{n) v^(n)/{i + (n) u + {n))\ converges to zero, cf. condition 
(|7.12[) . thus the right hand side in (|B. 16|) converges to zero for n — > oo. 

Finally, it remains to consider the term arising from the remainder in the expansion of q^, see 
(|B.13p . inserted into (|B.12|) . By assumption q'^ is bounded (by a constant q). Therefore we obtain 
an upper bound on the respective term by 



q\m\\L-- 



p(n) 
k=l 



U n {x) 



fc=l 

S 77 IU 112 j+M l 



vet 



L 1 ■ 



Here we have employed the Poincare inequality in L 1 with optimal Poincare constant given by 
diam(D fc , n )/2 [J. 

A combination of these estimates yields an upper bound to (|B.8[) by 
(ESI < C$ ^\\ Pj -z^(Q n )\\ L i + \\u-U^ L i+S + (n) 2 \\WU^ L i+5l(n)+5 + (n) + R(n)) , (B.17) 
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where the term R(n) is given by the right hand side of (|B.16[) and converges to zero for n — > oo. 
The constant C$ < oo is a suitable deterministic constant independent of n G N which depends 
on $ G (H 2s (D)) m via the norm in H S {D) of the components of $. 

(ii.2) Next we consider the mixed terms (|B.9[) . Analogously to part (ii.l) we apply the triangle 
inequality and obtain 



p 3 {x)qj l {u{x))(f)j{x)(j) i {x)dx - / z^(e n ){x)q jl (U n (x))(l} j {x)(j) i (x)dx 



+ E 



(n) Qfc,n 



q ji {U n {x))(j) j {x) c/)i(x)dx 



P( n ) gjfe," 



a «E77r^^ n (^")(^,i^,jB(^,i^,„>j 



^-^ l(k, n) 
fe=i v ' ; 

As in (ii.l) we obtain for the first term in this right hand side an upper bound by 

?||^IU-||^||£-|b i -^(e n )|| LI + L||^||£-||0 j ||L-||«-cr*|| £ i. 

Also the second term is treated as in (ii.l), i.e., applying the Mean Value Theorem and estimat- 
ing the resulting terms accordingly. In particular the remainder term is estimated completely 
analogously. Therefore, the only term we are left to estimate is 
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+ q\D k , n \ 



1 



D k J 



1 - 



4>i(x) 4>j(x) dx 



D k . 
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\D k J 



(j>i (x) dx 



D k 



1 



D, 



l(k,n)\D k>n \J \\D k 



4>i(x) dx 



\D k J 



D k ,n 

(x) dx 



4>j (x) dx 



First of all, using Young's inequality we obtain for the second term the estimate 



(IB~T9l) < 



1 



(B.18) 



(B.19) 



(B.20) 



which converges to zero for n — > oo. 

We next estimate the term (|B.18|) . Firstly, we note that as in part (a) we find using Poincare's 
inequality an upper bound to the term 



ID* 



1 



D 



k,n J D k 



i(x) - <pj(x)) 2 dx 



1 



D 



k,n J D k 



<pi(x) — 4>j(x) dx 



(B.21) 



and the upper bound is proportional to <5+(n) 2 . Next, expanding the two squared terms in (jB.21[) 
we find using the reverse triangle inequality that the term (|B.21I) is an upper bound to 



\D k , n \ 



D 



k,n J D 



(j)i(x) 2 dx + 



D 



k,n J D 



4>j{x) 2 dx 



\D k>n \ 



Di, 



-(77—/ (t>i(x)dx) ~ (-J— I 0j{x)dx) 

^J->k,n JD kt „ ' V Dfc,ra Jc t .„ J 

(j>i(x)^j(x)dx- (jj— / <l>i(x)dx}(— — / <f>j(x)dx 



Thus also this term possesses an upper bound which is proportional to S + (n) 2 . For n — > 00 the 
upper bound converges to zero. As for S + (n) — > also the term spanning the first and second 
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line converges to zero which was established in (ii.l), necessarily also the term in the third line 
converges to zero. Therefore we infer that the term (|B.18|) converges to zero proportional to 
5 + {nf. 

Now, a combination of these estimates yields analogously to ()B.17|) in (ii.l) that 

dBJl) < C*(\\ Pl - z™{e n )\\ L 2 + \\u-U n \\ L 2+5{nf\\VU n \\ L i+5 2 {n)+5{n) + R{n)y (B.22) 

Here R{n) is a term converging to zero for n — > oo arising from (|B.20[) and it is of the same type 
as the term R(n) in (ii.l). The deterministic constant C$ is independent of n £ N and depends 
on 3> via the norm in H S (D) of the components of $. 

(ii.3) A combination of the final results (|B.17[) and (|B.22|) in (ii.l) and (ii.2) yields that there 
exists a constant C$ < oo such that for almost all t it holds that 

$, G{u{t),p{t)) $> £ - a„($, G n (U?, Gl l ) <&> £ 

m 

< C*(^|b;(t)-2?(6niU= + Ki)-tf n (^ 



Here we have also employed the continuous embedding of L 2 (D) L 1 (D). We next square 
both sides of this inequality and integrate over (0,T). Afterwards we take the square root of the 
integral terms and further take the expectation of the resulting inequality. Finally, appropriate 
applications of Jensen's inequality yields that 

' $,G(u(t),p(t))$) s -a n ($,G n (U?,Q?)$) £ \dt (B.23) 

rn 

C*, T (^(n)+ ( 5H + i?(n)+E"[|| M -^ l || i2((o , T) , L2) +^|b l -zr(0 n )ll^((o,T),^) 

i=i 



o 



< 



for an appropriate constant CV,3> < oo. Note that in order to arrive at the estimate (|B.23|) we have 
further employed that the random term || VC/ rl ||x 2 ((o,r),i 2 ) can be estimated by a deterministic 
bound independent of n g N due to Proposition lB.il (b). Finally, due to the law of large numbers, 
i.e., Theorem l7.ll the sequence of PDMPs converges to the deterministic limit in the mean. Hence 
the expectation in the right hand side in (|B.23[) converges to zero for n — > oo. Furthermore, $+(n) 
converges to zero by assumption (|7.8[) . as does the term R(n). Thus, overall the right hand side 
in (|B.23I) converges to zero. The convergence of the quadratic variation is proved which completes 
the proof of Theorem 17.11 
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